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ABSTRACT 

This report was prepared to describe the results of 
an experimental study of free turbulent mixing in axial 
pressure gradients. Both the pressure field and the velocity 
field were extracted using pitot-static probes and both a 
single normal wire and an X-array hot-wire anemometers. 

The resulting hot-wire data shown in Figures 1-80 present 
an excellent visual representation of the physical aspects 
of turbulent flow. 

An analysis of the results of this experiment show 
conclusively that the Reynolds Stresses in the flow are 
related by some function to the turbulent kinetic energy; 
however, no specific function was found that would produce 
the relationship observed in this study. 
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NOMENCLATURE 



A 

T 


Boussinesq apparent viscosity 


A 


X-Array hot-wire lead (u+v) 


B 


X-Array hot-wire lead (u-v) 


C 


Constant in hot-wire calculations 


c i 


Constant in apparent viscosity relation 


c 

T 


Local turbulent shear coefficient 


C uv 


Correlation coefficient 


d 


Wake generator diameter 


e 


Linear output from hot-wire 


H 


Shape factor 


i J 


Indicates, i or j direction (i=l,2,... j = l,2,. 


k 


Turbulent kinetic energy 


S, 


Prandtl mixing length 


P 


Pressure 


q 


Dynamic pressure 


r 


General flow parameter 


S(x) 


Apparent viscosity constant 


T 


Characteristic time .* 


t 


Time 


x,y,z 


Spatial Coordinates 
x - Axial 
y - Lateral 


u, v,w 


Corresponding velocities 


V* 


Characteristic velocity 
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Constant in hot-wire calculations 
Displacement thickness 
Energy thickness 
Mixing length parameter = k™^ 

Momentum thickness 
Laminar viscosity 
Apparent viscosity 
Laminar kinetic viscosity 
Apparent kinetic viscosity 
Density 

Reynold's Stress 
Differential operator DEL 
Integral from a to b with respect to y 

Partial derivative with respect to r 

Ordinary derivative with respect to r 
Indicates absolute value 
Subscripts 

A Parameter of lead A on hot-wire 

B Parameter of lead B on hot-wire 

e Local external flow conditions 

x,y Indicates partial differentiation with respect to x 3 y 

t Shear related quantity 

Superscripts 

Indicates a time averaged quantity 
• Designates a fluctuating quantity 



8 

6 

6 3 

e 

e 

y 

v t 

V 

v t 

p 

T 

V 

b 

f dy 
a 

8u 

9r 

du 

dr 

I I 



6 



I . INTRODUCTION 



In most real fluid flows the transition from laminar to 
turbulent flow carries the fluid mechanist from a situation 
where he can safely predict the values of the flow-field 
parameters from the solution of the Navier-Stokes Equations 
to one where predictions are obtainable only on the average, 
at best. Yet most practical applications require at least 
some consideration of turbulence and how it can be dealt 
with. In the past the engineer has assumed various plausible 
hypotheses in order to introduce a measure of order and 
predictability to turbulent flow. It is the intention of 
this study, then, to attempt to experimentally verify certain 
of these hypotheses. Before describing the currently 
accepted assumptions regarding turbulent flow, a general 
review of their history will be presented to provide the 
necessary background information. 

A turbulent flow is defined as one in which an irregular 
fluctuation or mixing motion is superimposed on the princi- 
pal fluid motion. Its onset is accompanied by a noticeable 
change in the law of resistance. In the case of laminar 
pipe flow, the longitudinal pressure gradient which main- 
tains the motion is proportional to the first power of the 
velocity. By contrast, in turbulent flow this pressure 
gradient becomes nearly proportional to the square of the 
mean flow velocity. In order to calculate flow field values 
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in a turbulent field, then, some logical hypothesis must 
be made to account for this change in the apparent viscosity 
of the flow. First, it is generally assumed that the fluid 
motion can be separated into a mean and a fluctuating part 
which are linearly related to the total fluid motion in 
such a way that a flow parameter (r) can be defined 

r = r + r' (1-1) 

These definitions can be substituted into the equations of 
motion to produce a new set of equations involving mean and 
fluctuating terms and their cross products. Assuming that 
the time average of the fluctuations is zero, the equations 
may be averaged to get the following momentum equation for 
a two-dimensional mean flow (see Appendix A for the complete 
development ) : 



u 



-3u 

9x 



— 3u 

v 9y 



- u 



du, 
e dx 



+ v, (x) 



3 2 u 

3y 2 



( A-22 ) 



where the apparent kinematic viscosity (v t (x)) is a function 
of the turbulent flow field, and is calculated by using some 
logical empirical hypothesis for approximating its value. 

Since only the mean values of the flow parameters are 
known, some means of relating these values to the Reynolds 
stresses produced by the fluctuating part of the flow is 
required. The first hypothesis for relating these quantities 
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was proposed by J. Boussinesq, 



t = -pu'v' = A (1-2) 
o y t 3y 

where is an apparent or eddy viscosity that is not a 
fluid property but depends on the structure of the flow 
field. 



A t = p v fc (x) 



This apparent viscosity is generally assumed to be propor- 
tional to some characteristic velocity (v’"’) times the mean 
velocity gradient: 



A 



T 



V 



* 5u 

9y 



(1-3) 



In 1925 Prandtl made an important advance in turbulent 
theory by introducing the concept of mixing length and the 
relation, 



T 



t 




du I du 
dy 1 dy 



u-*n 



v t (x) = l 2 




where 



l = mixing length 



Z 




characteristic velocity 
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This relation has produced relatively accurate values for 
some limited cases. However, this definition of the kine- 
matic viscosity implies that it should vanish where the 
velocity gradient is zero, v;hich is contradictory to experi- 
mental evidence. Several subsequent hypotheses were proposed 
by Taylor, Von Karman, and Prandtl himself, but they were 
all limited in application to free turbulent flow. 

Later, in an effort to develop a universally valid 
model for the function y^, Prandtl introduced an energy- 
based hypothesis which was subsequently improved on by 
Emmons ( 195 * 0 , and Glushko ( 1965 ), 

y t - p k 1/2 a (i-5) 

k = (u 2 + v 2 )/2 

Here k was found to be a better characteristic velocity than 
in the earlier mixing length hypothesis, but the length 
scales were difficult to ascertain. More modern theories, 
referred to as two-equation viscosity models, use sophisti- 
cated analytical methods to determine this length scale. 

The apparent viscosity is first assumed to be related to the 
kinetic energy. Then some hypothesis of the form 
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is made to make it possible to calculate the characteristic 
length. One of the most accepted of these relations was 
introduced by Chou in 19^5: 




( 1 - 6 ) 



2 

M t = C 1 p k /e 



Thus 




( 1 - 7 ) 



It was this type of assumption which initiated the pro- 
ject reported here. Using both pitot-static and hot-wire 
anemometer sensors, an attempt was made to validate the 
hypothesis that shear was predominantly a function of the 
turbulent kinetic energy. 
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II . EQUIPMENT 



A . WIND TUNNEL 

1 . General Description 

The tunnel consists essentially of three major 
components: the inlet, the test section, and the blower as 
shown in Figures 81-83. The inlet was designed to contract 
the slow moving inlet air to a high-velocity parallel flow 
at the entrance of the test section. At this point, a 
Plexiglas cylinder was introduced into the flow to generate 
a wake velocity defect which was then subjected to a constant 
adverse pressure gradient. After leaving the test section, 
the flow w as maintained at a constant area through a length 
of tubing, then diffused and directed into the centrifugal 
blower and exhausted vertically. 

2 . The Inlet 

The inlet section of the tunnel was composed of the 
inlet bellmouth, the plenum chamber, and the nozzle. Incoming 
air was first contracted into a two-foot square plenum sec- 
tion. This contraction was accomplished by the inlet bell- 
mouth constructed from eight-inch diameter aluminum pipe 
halves. The corner seams were filled with a sanded epoxy 
fillet to effect a smooth transition of the air into parallel 
flow. The plenum section consisted of four two-foot square 
by one-foot long steel sections bolted together. Cork gas- 
kets were inserted between the sections to prevent the influx 
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of air at the joints and to allow for undisturbed air flow 
over these junctions. A four-inch thick section of quarter- 
inch honeycomb was installed to straighten the air flow and 
prevent swirling. Next, two copper screens were inserted 
to reduce the ambient turbulence level in the flow. Finally, 
the air passed from the plenum chamber to the test section 
through a two-foot long steel nozzle with an exit area of 
forty-nine square inches (seven inches square), and produced 
a contraction ratio of 11.76. Static pressure ports were 
located on the right side of the nozzle at the inlet and 
exit for determination of the exit dynamic pressure. 

3 . The Test Section 

The test section consisted of two parallel, vertical 
side walls separated by flexible and adjustable top and 
bottom- walls . The flexible walls were constructed from 
one-eighth inch Plexiglas fastened to the nozzle with piano 
hinges. As described in Reference 1, these walls could be 
varied to produce a variety of different pressure gradients, 
but for this experiment they were set to create a constant 
adverse pressure gradient followed by a constant favorable 
pressure gradient. 

One vertical wall was provided with thirteen slots 
for inserting data sensors. These slots were nineteen 
inches high, three-eighths of an inch wide, and spaced five 
inches apart. A filler plug for each slot was carefully 
hand milled, and polished to prevent disturbing the flow 
along the wall when the slots were not in use. A thin 
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aluminum strip was wedged over each filler plug to ensure a 
tight fit. When these plugs were removed to insert probes, 
a sliding teflon seal was used to prevent an influx of air 
through the slot (see Figure 8*1). 

The other vertical wall of the test section had three 
rectangular access holes along the axial centerline for 
adjustments and measurements. In addition, as one can see 
in Fig. 85, thirty-three static pressure ports were inserted 
along the axial centerline at one-and-one-half inch inter- 
vals. A manometer bank, shown in Fig. 86, connected to these 
static pressure ports provided a visual indication of the 
uniformity of the pressure distribution, as well as indica- 
ting the presence of any regions of separation. As an 
additional indication of possible flow separation, several 
rows of tufts were attached to the walls in critical areas. 

Having produced the desired environment for the 
type of air flow that this study was initiated to investi- 
gate, a turbulent wake velocity defect was introduced into 
the tunnel by placing a one-inch diameter cylinder in the 
inlet of the test section. In order to ensure a turbulent 
boundary layer at separation from the cylinder, a thin strip 
of abrasive, shown in Fig. 87» was attached to the stagnation 
line of the cylinder. 

4 . The Blower 

After exiting the test section, the air passed 
through two seven-inch square, four-foot long steel sections 
into a diffusor and then into the blower. A centrifugal 
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blower operating at three thousand revolutions per minute 
was utilized to pull air through the tunnel. The air was 
then exhausted vertically through sections of steel pipe 
which contained two dampers for the control of the total 
mass flow through the tunnel. 

B. DATA SENSORS 

1 . General Description 

The data for this study were taken with the three 
different types of probes shown in Fig. 88: A pitot static 

probe, a single normal hot wire probe, and an x-wire hot 
wire probe. A detachable traversing mechanism was used to 
position the probes vertically, and the values were recorded 
by hand or on an x-y plotter. 

2 . The Traverse 

This apparatus was equipped with two constant speed 
motors to drive the mechanism across the nineteen-inch 
vertical span of the tunnel. The probe was clamped to a 
movable support on the traverse, and the vertical position 
measured by indexing a reference mark on this mechanism to a 
scale fixed to the frame of the apparatus. This scale was 
readable to one one-hundreth of an inch. A rack-and-pinion 
gearing system drove a potentiometer that was also attached 
to the movable support to produce an electrical position 
signal for the abscissa of the x-y plotter. 
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3 • The Pitot-Static Sensor 



A Prandtl-type pitot-static pressure probe Model 
# PBA-12-F-11-KL, built by United Sensors and measuring 
8l hundreths of an inch from stem to tip, was used to obtain 
dynamic and local static pressures at one-tenth of an inch 
intervals through the wake. These pressures were read on 
two variable angle, water-f illed U-tube manometers referenced 
to atmospheric pressure. This reference was chosen to re- 
duce the fluctuations in the water levels. As described in 
Reference 1, taking P-total minus P-static directly creates 
very large fluctuations that make accurate data taking diffi- 
cult. By separating the two readings, the fluctuations in 
the manometer were significantly reduced, but not 
eliminated . 

A third manometer was attached to the plenum chamber 
of the tunnel, and referenced to the last static pressure 
port on the test section. As one can see in Fig. 89 , all of 
these manometers were equipped with portable plastic sights 
which restricted one's line of sight to a plane perpendicu- 
lar to the manometer, thus precluding possible parallax and 
assuring consistent readings. 

4 . The Normal Wire 

This sensor, shown in Fig. 90, was a single-wire 
anemometer probe, Model § 1150 built by Thermo-Systems Inc. 
(TSI), measuring 4 inches from stem to tip. Again the probe 
was fixed to the traverse for data collection. The lead 
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from the sensor was connected to a 105^A-30 anemometer 
module with a 1051-2 monitor and power supply, also built 
by Thermo-Systems Inc. (TSI) (See Fig. 91 & 92), which 
converted the voltage changes to readable data. Next these 
data were fed into a TSI 1015C correlator that was used to 
separate the alternating and D.C. signals (see Appendix B). 
The output of the correlator was connected to an analog 
data processing computer which indicated the mean square 
average output of the correlator. Finally, the output of 
this apparatus was connected to the ordinate of the X-Y 
plotter, recalling that a potentiometer on the traverse was 
connected to the abscissa. 

5 . The X-Wire Anemometer Probe 

The last sensor, a Model ff 1055 built by TSI, con- 
sisted of two hot wires oriented at ^5° to the axis of the 
probe such that the output signals of the wires corresponded 
to u + v and u - v respectively (See Appendix B) . These 
data were processed in the same manner as the normal wire 
data except for the use of several additional algebraic 
functions of the 1015C correlator. 

j 
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III. EXPERIMENTAL PROCEDURE 
A. GENERAL DESCRIPTION 

The previous study described in Reference 1 had con- 
figured the test section to produce a constant, adverse- 
favorable pressure gradient. Since this was the environment 
required for this study, no changes were made in the test 
section except to ensure a tight seal between the vertical 
and the movable walls by the addition of strips of clay 
along the lines of contact. 

The actual collection of data was carried out in three 
distinct steps. First, the pitot-static system was used to 
collect the total and static pressure data that were intro- 
duced into the computer program, described in Appendix E, 
and used to calculate the flow parameters. Next, a normal- 
wire anemometer was inserted in the test section to gather 
information on both the mean and fluctuating axial velocities. 
These hot-wire values were converted to velocities by equa- 
ting the maximum voltage at each station to the maximum value 
of the velocity calculated from the pitot-static data. 

Finally, an X-array hot-wire anemometer probe was used to 
get both mean and fluctuating data for velocities in both 
the axial and the lateral directions, as well as Reynolds 
stresses. Throughout the experiment a constant wind tunnel 
speed was maintained by ensuring that the total pressure of 
the Plenum chamber, minus the static pressure from the test 
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section, was no less than 9.90 nor greater than 10.00 
inches of water on a thirty-degree inclined manometer. 

B. PITOT-STATIC DATA 

In this phase of the experiment all data were recorded 
(to an estimated accuracy of one one-hundreth of an inch) 
from water-filled U-Tube manometers connected to the pitot- 
static probe. At each of eleven stations downstream from 
the cylinder the static and total pressures were recorded 
at intervals through the wake of one-tenth of an inch. At 
each station the probe was moved near the upper limit of 
the traverse and situated in the free stream, then the tra- 
verse was marched across the flow, recording the values of 
total and static pressure at each increment. The initial 
determination of the freestream was accomplished by observing 
that the manometer becomes almost totally insensitive to 
changes in the vertical direction in the undisturbed flow. 

A change in total pressure of less than three one-hundredths 
of an inch for three consecutive positions was taken to be 
the edge of the wake defect. 

As the traverse approached the center of the wake the 
intensity of the turbulence became so large that accurate 
readings were difficult to obtain due to fluctuations in the 
water level. However, by observing the water levels for a 
long period of time one could determine an average height 
that is estimated to be accurate to three one-hundredths of 
one inch. The same criterion was used at the bottom of the 
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flow channel, as at the top, to determine when the 
freestream had been reached. 

C. AXIAL FLOW ANEMOMETER DATA 

The probe used in this phase was a single normal hot- 
wire anemometer designed to produce voltage changes corres- 
ponding to velocity fluctuations in the flow. VJhen the 
tungsten wire was exposed to the air in the wind tunnel it 
was cooled by the flow over it, thus requiring additional 
voltage to maintain a constant temperature in the wire. It 
is apparent, then, that the velocity in the tunnel is pro- 
portional to the voltage required to maintain a constant 
temperature in the hot-wire. The probe is placed in a 
bridge circuit in the TSI module, and its output linearized 
to produce a linear function of velocity of the form: 

e = c (u + u') (III-l) 

where , 

u = u + u* 

u = mean flow velocity 
u* = fluctuating velocity 

and 

T 

/ u' dT = 0 ( III-2 ) 

0 

Therefore if the DC output on the correlator is selected, 
the output would be 
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e = c (u + u ' ) 



However, since this signal was a strong function of time, 
the data were averaged using a low pass filter with a time 
constant (T) equal to 0.6 secs. The final output of the 
averager was simply 



e = cu (III-3) 

and this quantity was plotted versus the traverse position. 

The second type of data for this phase was taken with 
the selector switch of the correlator on AC, thus producing 
an output corresponding to the fluctuating portion of the 
flow : 



e* = cu’ (III-4) 

This signal was also a strong function of time but, because 
of the nature of the term u’, as expressed in (III-2), the 
signal could not be directly averaged. Therefore, in order 
to provide a finite measure, u' was first squared and then 
averaged and plotted versus the traverse position, since 



T 

is / (u* ) 2 dx t 0 

1 0 
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Once a plot of e and e' had been obtained, the velocity 
could be calculated by solving equations (III-3) and (III-4) 
to get 



u = 1/C e 




where from equation (B-5) in Appendix B 

6 

C h — (III-5) 

u 

e 

In order to produce accurate values for u and u' in the 
laboratory environment the equipment had to be carefully 
adjusted for each individual flow. Specific instructions 
for this procedure can be found in the TSI operational 
manual but, in general, the first step was to adjust the 
control resistance to balance the bridge circuit.. This 
adjustment was accomplished by setting the control resistor 
to a value that would produce a specified voltage change on 
the monitor. At this point the system was ready to run and 
the next step was to select the linearized operational mode 
on the monitor, position the probe in the free stream, and 
maximize the frequency response of the wire by setting the 
inductance to the minimum value that would not cause the 
signal to oscillate. The stability of the system was 



22 



monitored during this step on a dual-beam oscilloscope to 
ensure no oscillations were present in the signal. When 
the probe had been adjusted to the flow a digital voltmeter 
was used to check the multiplier-averager for correct opera- 
tion. Then the x-y plotter was set to record the signal on 
an appropriate scale. 

Now that the system had been checked for correct operation, 
it was positioned as close to the upper wall as practicable 
and the pen was dropped on the plotter. The sensor was then 
driven across the wake to a point near the bottom wall, and 
the pen lifted. After changing the paper in the plotter, 
the same procedure was used to make a traverse back to the 
top of the test section. A total of four traverses were 
made at each station to record both the mean flow and the 
mean square of the fluctuating flow. Each of these quan- 
tities was recorded at two different scales in order to 
facilitate translation into digital data for the IBM-360. 

D. TWO-DIMENSIONAL ANEMOMETER DATA 

The basic operation and accuracy checks of the x-array 
hot wire are very similar to those of the single normal 
wire. However, some additional complications are created 
by having two separate wires on the probe. The terms are 
now defined. 



e = e + e ’ 
u = u + u' 
v = v + v* 
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As explained in Appendix B, the outputs of the two wires 
(A & B) were 



e A “ C A <u + v) 
e B ” °B (u ' v) 

Thus, for a given flow, the values of C. and C- cannot be 

A ID 

determined easily. In order to find these constants one 
must have 



In this study the author used a basic physical character- 
istic of pipe flow to set these constants equal, i.e., that 
the lateral velocity must be zero at the centerline. Recall 
that, at y = 0 



'A 



■ °A U ' 



C B U = 



'B 



Thus the calue of (e.-e D ) was set equal to zero so the 
constant (C) could be defined as in (III-5). Note that this 
procedure is equivalent to setting V = 0 at the centerline 
since. 



e A - e B = C(u + v) - C(u - v) = 2Cv 



2k 



e A + e B = + v ) + C(u - v) = 2Cu 



The actual position of the centerline was determined from 
the pitot-static data taken in phase 1. 

Again, the wires were adjusted to the individual flow at 
each station and the data were taken in the same manner as 
described in the previous section. However, several more 
traverses were required to get all of the data available 
to the x-array hot wire. The settings on the correlator, 
and the corresponding flow parameters, are listed below. 

All barred quantities were taken with the correlator on 
DC, and averaged before plotting. The primed quantities 
were taken with the correlator on AC, and squared, averaged, 
and then plotted. 



SETTING 



(A+B) 

(A^B) 

A 



B 



RESULT 



u 



V 



u+v 



u-v 



SETTING 



(A+B) ‘ 2 



(A-B) ,2 




RESULT 




(u+v) • 2 



(u-v) ' 2 
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IV. DISCUSSION OF RESULTS 



Since this study was intended to carry out an investi- 
gation of free turbulent mixing in a constant pressure 
gradient, it would be appropriate to discuss how close the 
actual experimental environment approached that desired for 
the study. It can be easily seen in Fig. 21 that a wake 
velocity defect was initially created by the cylinder as 
desired, and the magnitude of u* in Fig. 4l clearly shows 
that the flow v/as, in fact, turbulent. However, it can 
also be seen that the adverse pressure gradient encourages 
the formation of a large boundary layer which, as it mixes 
downstream, interacts with the wake defect, tending to 
raise some doubt that free turbulent mixing is still present 
downstream of Station 6. It is especially significant to 
note that the turbulence near the wall builds up from a 
negligible value at Station 1 (see Fig. ^1-50) to a value 
greater than that induced into the wake by the cylinder 
midway down the test section. As expected, the wake defect 
eventually begins to mix out 'in the favorable pressure 
gradient. 

By comparing the stress calculated from the pitot-static 
data with that from the hot-wire anemometer, one can easily 
see that in both cases the stress is zero on the centerline, 
then builds up to a maximum as the velocity gradient nears 
the region of its maximum. At the later stations, both 
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types of data show a shear stress sign reversal, and its 
magnitude begins to increase as the traverse approaches the 
edge of the test section. This trend in the data was prob- 
ably reflecting the increasing interaction with the boundary 
layer as the air flowed downstream. 

Due to the greater capabilities of the hot-wire probe, 
it was used to take data very near the wall, and the velocity 
can be seen to fall rapidly in the vicinity of the boundary 
layer. The hot-wire data (Figures 1-80) were then used to 
prepare Table I to display the value of the correlation 
coefficient (C uy ) and n in the turbulent flow field. Note 
that the correlation coefficient C uy approaches zero at the 
centerline, and reaches its maximum where the shear is a 
maximum as predicted by Schlicting (See Fig. 18.^1 Ref. 2). 

Considering the relation between kinetic energy and the 
shear stress, it is easy to see in Table II that both param- 
eters reach a maximum at nearly the same point in the flow. 
Comparing the anemometer with the pitot-static data in 
Appendix F, the maximum shear there also occurs at the 
lateral position corresponding to the kinetic energy maxi- 
mum. Thus, it is clear that at least some significant func- 
tional relationship exists between them. The ideal relation- 
ship would be that postulated by the theory of eddy viscosity, 
that the two parameters differ only by a constant that is 
a function of x. The slight displacement of the two maxi- 
mums in this case would imply that the relationship between 



/ 
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k and is, in reality, a more complex function. As an 
additional attempt to clarify the relation between kinetic 
energy and shear stress, equation (7) was approximated by 
taking the characteristic velocity, as the square root of 
the kinetic energy, and 



s = c 1 a = s(x) (iv-i) 

As a result, the relation for Reynolds stress becomes, 



3u 

9y 



k = y*. ^ = S(x) k 



. 1/2 



3u 

3y 



( IV- 2 ) 



The values for equation (IV-2) were calculated and plotted 
versus traverse position for several stations from the data 
in Table II> which was derived from Figures 21-80. The 
results of these plots are shown in Figures 93-97. Once 
again, there is a definite similarity between the results 
of the approximating equation, (7), and the calculation of 

_________ S 

u'v' from the hot wire data. In this case the shapes of 
the curves were exceptionally similar, but the maximums 
were still displaced by a small amount, preventing direct 
correlation by multiplication by a constant. The author 
feels that a more comprehensive analysis of the data in 
this report is warranted by the obvious visual relationship 
of the curves. An in-depth analysis might well uncover 
a significantly better model for calculating turbulent 
flow fields than are currently in use. 
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TABLE I 



Calculation of The Correlation 
Coefficient and Eta 



/'•' .... o: : - n== tief+ej* ) 

Al/oZo x' 2 - 1 \!/^ TvT ^ 5 I t ^ ^ \< * i / /-■. 



^ fe+e')' x 1( 


&-er ^ (e A teJ x -Ke^-a^r 


Station 1 


Pnsi -hi nn _ 




a. 


10.^ 


-.6529 


\ 

1.310 


11.0 


-.4881 


1.341 


11.5 


-.4353 


1.494 


12.0 


-.2092 


1.745 


12.5 


-.0319 


1.812 


13.0 


.1501 


1.782 


13.5 


.2857 


1.489 


lit-.o 


.3112 


1.292 


14.5 


.1631 


1.243 




1 

Station 2 




9.5 


-.4045 


.9032 


10.0 


-.4170 


1.250 


10.5 


-.2029 


1.012 


11.0 


-.2353 


.9706 


n.5 


-.2470 


1.058 


12.0 


-.3108 


1.243 
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Posi tion 


C 


n 








12.5 


-.0381 


1.219 


13.0 


.1172 


1.249 


13.5 


.1880 


1.188 


llj-.O 


.25 76 


1.343 


1^.5 


.0851 


1.277 


15.0 


.0683 


1.153 


15.5 


-.2887 


.875- 


Station 3 


9.0 


- . 4544 


.878 


10.0 


-.3404 


.922 


11.0 


-.3861 


1.122 


12.0 


-.2401 


1.227 


13.0 


.0915 


1.281 


14.0 


.1889 


1.211 


15.0 


.1097 


1.071 


Station 4 


8.0 


-.6876 


1.108 


9.0 


-.5699 


1.281 


10.0 


-.4562 


1.361 


11.0 


-.3599 


1.426 


12.0 


-.1709 


1.472 


13.0 


.0983 


1.474 


14.0 


.2615 


1.485 
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4 (cont.) 




Position 


C 


_ A 


15.0 


.2815 


V 

1.314 


16.0 


.2400 


1.194 


17.0 


-.0423 


1.036 




Station 5 




8.0 


-.2750 


1.083 


9.0 


-.3340 


1.176 


10.0 


-.3470 


1.317 


11.0 


-.3040 


1.341 


12.0 


-.1508 


1.388 


13.0 


.0582 


1.427 


14.0 


.1907 


1.390 


15.0 


.2584 


1.317 


16.0 


.1848 


1.161 


17.0 


-.0416 


1.105 




Station 6 




7.0 


-.0154 


1 . 145 


8.0 


-.2951 


1.125 


9.0 


-.2809 


1.165 


10.0 


-.2711 


1.269 


11.0 


-.4096 


1.119 


12.0 


-.0900 


1.254 


13.0 


.0295 


1.374 
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Station 8 


Position 


c 


V ---- - 


7.0 


.1517 


\ 

1.144 


8.0 


-.1027 


1.010 


9.0 


-.3392 


1.138 


10.0 


-.7895 


1.095 


11.0 


-.7298 


1.194 


12.0 


-.3909 


1.205 


13.0 


.1279 


1.165 


14.0 


.4510 


1.232 


15.0 


.2654 


1.200 


16.0 


.0634 


1.151 


17.0 


-.4380 


1.185 


Station 10 


9.0 


.0454 


1.095 


10.0 


-.7389 


1.130 


11.0 


-.8439 


1.079 


12.0 


-.3313 


I.080 


13.0 


.2910 


1.091 


lit. 0 


.4818 


1.105 


15.0 


-.0229 


1.145 


16.0 


-.0116 


1.020 
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Station 9 










Position 


C 


> 


9.0 


-.3614 


\ 

1.123 


10.0 


-.6485 


1.063 


11.0 


-.7565 


1.130 


12.0 


-.3710 


1.091 


13.0 


.1835 


1.143 


14.0 


.5674 


1.146 


15.0 


.3130 


1.140 
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TABLE II 



Calculation of Kinetic Energy 
And Shear Stress Parameters 




Station 1 


Position 


k 


u'v' 


12.0 


993.9 


36,81 


13.0 


969.0 


70.99 


13.2 


920.2 


84.14 


13.4 


O70.3 


94.66 


13.6 


780.9 


115.69 


13.0 


697.3 


89.40 


14.0 


912.7 


78.88 


14.2 


331.3 


47.33 


14.4 


239.3 


34.18 


14.6 


199.1 


. 13.15 




Station 2 




12.0 


939.0 


-83.0 


12.5 


945.0 


10.4 
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Station 2 ( cont . ) 


Position 


k 


u * v ' 


12.8 


534.0 


20.74 


13.0 


531.4 


31.11 


13.2 


521.0 


36.29 


13.4 


505.5 


31.11 


13.8 


J4 66. 6 


20.74 


13.8 


414.8 


7.78 


14.0 


362.9 


5.18 


14.2 


3H.1 


2.59 


14.4 


267.0 


5.18 


14.8 


223.0 


5.18 


14.8 


189.3 


7.78 


15.0 


153.0 


1 5.18 


Station 3 


12.4 


378.5 


-22.08 


12.6 


376.9 


-11.04 


12.8 


376.9 


3.15 


13.0 


380.1 


15.77 


13.2 


381.7 


28.39 


13.4 


378.5 


28.39 


13.6 


364.3 


29.97 


13.8 


343.8 


29.97 
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Station 3 (cont 


.) 


Position 


k 


u ' v 1 


14.0 


320.2 


29.97 


14.2 


296.5 


23. 66 


14.4 


265.0 


22.08 


14.6 


236.6 


17.35 


14.8 


205.0 


12.62 


15.0 


176.6 


9.46 


15.2 


135.6 


6.31 


15.4 


119.9 


3.15 


Station 4 


13.0 


235.9 


11.53 


13.2 


237.5 


16. 14 


13.4 


238.2 


21.52 


13.6 


237.5 


25.36 


13.8 


235.1 


27.66 


14.0 


232.8 


29.97 


14.2 


226.7 


32.28 


14 • 4 


219.0 


30.74 


14.8 


209.8 


25.36 


14.8 


199.0 


23.05 


15.0 


186.0 


19.98 


15.2 


172.9 


16.91 



37 



— — — 

Station 4 (cont.) 


Position 


k 


u 1 V' 


15.4 


156.8 


16.91 


15.6 


139.9 


12.30 


15.8 


139.1 


10.76 


16.0 


101.4 


8. 45 


16.2 


86.1 


7.69 


16.4 


72.2 


6.15 


16.6 


61.5 


3.84 


16.8 


52.3 


3.84 


Station 5 


12.0 


217.5 


-14.32 


12.2 


214.8 


- 8.05 


12.4 


213.0 


- 3.58 


12.6 


213.0 


0.90 


12.8 


215.7 


4.48 


13.0 


217.5 


7.16 


13.2 


221.0 


11.63 


13.4 


224.6 


14.32 


13.6 


228.2 


17.00 


13.8 


230.0 


19.69 


14.0 


232.7 


21.43 


14.2 


231.3 


24.16 


14.4 


• , 229.1 


25.95 
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Station 5 (cent 


.) 


Position 


k 


U' V' 


II4.. 6 


226.4 


25.95 


14.8 


221.9 


25.95 


15.0 


214.8 


26.85 


15.2 


208.3 


26. 85 


15.4- 


197.3 


26 . 85 


15.6 


187.9 


24.16 


15.8 


176.3 


21.48 


16.0 


166.5 


17.90 


16.2 


152.1 


14.32 


16.4 


137.8 


11.63 


16.6 


125.3 


8.95 


16.8 


107.4 


1.79 


17.0 


93.97 


- 1.79 


17.2 


85.0 


- 2.69 


17.4 


78.8 


- 3.58 




Station 6 


- 


13.0 


191.1 


2.71 


13.2 


197.5 


6.34 


13.4 


202.0 


10.87 


13.6 


205.6 


13.59 


13.8 


208.4 


18.12 


14.0 


211.1 


20.84 
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Station 6 ( cont. ) 


Position 


k 


u 1 v 1 


14.2 


212.0 


23.55 


14.4 


212.0 


24.46 


14.6 


211.1 


24*48 


14.0 


208.4 


23.55 


15.0 


204.7 


22,65 


15.2 


199.3 


21.74 


15.4 


193.9 


22.65 


15.6 


188.4 


19.02 


15.8 


181.2 


16.31 


16.0 


174.8 


13.59 


16.2 


164.0 


10.87 


16.4 


153.1 


7.25 


16.6 


144.9 


5.44 


16.6 


135.0 


3.62 


17.0 


129.5 


2.72 


17.2 


121.4 


0.0 


17.4 


116.9 

1 


- 3.62 


Station 7 


13.0 


169.3 


10.37 


13.2 


173.9 


14.49 


13.4 


178.4 


17.21 


13.6 


182.9 




18.11 
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Station 7 (cont 


.) 


Position 


k 


U ’ V 1 


13.8 


188.4 


18.11 


34.0 


192.9 


18.11 


14.2 


195.6 


18.11 


14.4 


198.3 


18.11 


14 • ^ 


198.3 


18.11 


14. s 


215.0 


18.11 


15.0 


193.8 


17.21 


1^.2 


190.2 


17.21 


15.4 


184.7 


15.^0 


15.6 


180.2 


14.49 


15.8 


173.9 


11.77 


16.0 


169.3 


9.06 


16.2 


163.9 


8.15 


16.4 


158.5 


7.25 


16.6 


154.0 


6.34 


16.8 


151.3 


4.53 


17.0 


148.5 


1.81 


17.2 


147.6 


- 3.62 


17.4 


149.4 


-10.87 


17.6 


154.9 


-16. 30 
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Station 8 


Position 


k 


u'v« 


15.0 


171.8 


11.34 


1^.2 


170.1 


9.60 


13.4 


166. 6 


9.60 


15.6 


164.0 


7.83 


15.8 


139.6 


3.23 


16.0 


136.1 


2.62 


16.2 


131.8 


- 1.73 


16.4 


148.3 


- 3.23 


16.6 


147.4 


- 9.60 


16.8 


130.0 


- 14.83 


17.0 


160.3 


- 18.32 


Station 9 


12.0 


136.4 


-13.22 


12.2 


133.1 


- 9.92 


12.4 


131.4 


- 3.79 


12.6 


131.4 


- 1.63 


12.8 


133.9 


2.48 


13.0 


138.8 


6.61 


13.2 


140.3 


10.74 


13.4 


143.8 


14.03 


13.6 


147.1 


17.36 


13.8 


131.2 


20.66 
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Station 9 (cont 


.) 


Position 


k 


u»v» 


14.0 


152.9 


21.49 


14.2 


152.1 


20.66 


14.4 


11(9.6 


19.83 


14.6 


147.1 


16.53 


14.8 


143.8 


14.88 


15.0 


142.1 


11.57 


Station 10 


12.0 


123.0 


-10.66 


12.2 


119.7 


- 6.56 


12.4 


119.7 


- 1.64 


12.6 


121.4 


2.46 


12.8 


124.6 


6.56 


13.0 


126.3 


9.84 


13.2 


132.0 


13.94 


13.4 


135.3 


16.40 


13.6 


138.6 


18.04 


13.8 


141.1 


19.68 


14*0 


140.2 


17.22 


l4o2 


138.6 


15.58 


14.4 


138.6 


12.30 


14.8 


137.8 


7.38 


14.8 


136.1 


3.28 
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Station. 1 Normal Wir© Moan Velocity 



ft/volt- S_££ 




o 

r— I 



oo o U\ 

(S4X oa ) H - ^T°°I e A uboh 



to 




G 




X2 




o 




d 




•H 




' — ' 


OJ 


d 




o 






r£j 


-P 




•H 


!ro 


CO 


O 


o 


M 


CL, 




G 




CO 




U 




o 




> 












EH 





co 



^5 



tation. 2 Normal Wire Mean Velocity 




Ile64 ft/volt- 




r n 



i 



l 



e. 






n 



n±t 









TTT- 



ir~ 



d-H 



P 



r~ 



m 






x= 






Trfr 



, rg qp 



46 



33 






i 



EEt 



as 



±t 



P 



4- 












~ r 1 



•“{tt 



-'hr 1 ' 



#r 



6334 






a 

rnx 



-f-Tffi 

t T!ir 



m 



?•. 



I 



tip 



t 



+f^ 



; S£ 



dfr 



44 



.:rr-1 



33 



o 

* 

v£) 

H 



53 



4xrf 

n±d 



^4- 



33 



'HfT- 






Hff 



Ttfl 

; rp‘ 

p 



P 



- j a 

Tpr n±r 



O 

a 

r-i 



o 

OJ* 



rtt£ o 
S4 h 



pit 



O 

o 



O CO 

rH _ i — I 

(64X OA ) 3 - £4T30T 9 A UB9K 



o 

1A 



co 

G> 

-d 

O 

a 

*H 



d 

o 

•H 

-P 

•H 

CO 

o 

cu 

a) 

00 

d 

G) 

> 

ctf 

d 



CA 



Kl 

cd 

o 

M 



46 



Station 3 Normal Wire Mean Velocity 
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FIGURE 5 

Station 5 Normal Wire Mean Veloci 
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Station. 8 Normal Wire Mean Velocity 



26.65 i't/volt-sec 
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FIGURE 9 

Station 9 Normal Wire Mean Velocity 
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FIGURE 10 

Station 10 Normal Wire Moan Velocity 



= 157.5 ftVvolt-sec 
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FIGURE 11 

Station 1 MS Normal Wire Fluctuating Velocity 



136,3 1't /volt-sec 
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FIGURE 12 

Station 2 MS Normal Wire Fluctuating Velocity 
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FIGURE 13 

Station 3 MS Normal Wire Fluctuating Velocity 
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FIGURE 14 

Station 4 MS Normal Wire Fluctuating Velocity 
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FIGURE 15 

Station 5 MS Normal Wire Fluctuating Velocity 
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629#7 ft /volt-sec 
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Travers© Position (inches) 

FIGURE 17 

Station 7 MS Normal Wire Fluctuating Velocity 
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Station 1 X-Array Mean Axial Velocity 



22. 77 ft/volt-sec 
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FIGURE 22 

Station 2 X-Array Mean Axial Velocity 
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Station 3 X-Array Mean Axial Velocity 




12.U0 ft/volt-sec 
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FIGURE 214- 

Station 4 X-Array Mean Axial Velocity 



13.38 ft/volt- 
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FIGURE 25 

Station 5 X-Array Mean Axial Velocity 



13. 46 ft/volt- 
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13*46 ft/volt- 
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FIGURE 2 7 

Station 7 X-Array Mean Axial Velocity 




13.21 ft/volt-sec 
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FIGURE 28 

Station 8 X-Array Mean Axial Veloci 



12.86 ft/volt-sec 
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FIGURE 29 

Station 9 X-Array Mean Axial Velocity 



12.81 ft/volt-sec 
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FIGURE 30 

Station 10 X-Array Mean Axial Velocity 



22.93 ft/volt-sec 
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FIGURE 31 

Station 1 X-Array Mean Lateral Velocity 



22.77 ft/volt-sec 
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FIGURE 32 

Station 2 X-Array Mean Lateral Velocity 



ft /volt 
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+ 1.5 



j3 = 12.40 ft/volt -sec 
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Position (inches) 



13.33 ft/volt-sec 



fffrr 




o 




* 


CO 




o 


nA 


xJ 




o 




£ 




•H 








a 


O 


O 


e 


•H 


C\i 


■P 


H 


•H 




OQ 




o 




Ph 




<d 




CQ 


O 


Pi 


« 


o 


O 


> 




cd 




Pi 




Eh 



O 

CO 



o 

vO 



(S4X0A) a - ^q.T00i0A TBJoq.B'i ubgh 



78 



FIC-URS 35 

Station 5 X-Array Mean Lateral Velocity 




13 . 4-6 ft/volt-sec 
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FIGURE 36 

Station 6 X-Array Mean Lateral Velocity 





13.ij-6 ft/volt 
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FIGURE 37 

Station 7 X-Array Mean Lateral Velocity 




13.21 ft/volt-sec 
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12.86 rt/volt-sec 
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FIGURE 39 

Station 9 X-Array Mean Lateral Velocity 




12.81 ft/volt-sec 
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FIGURE Lj.0 

Station 10 X-Array Moan Lateral Velocity 
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FIGURE J+l 

Station 1 X- Array MS Fluctuating Axial Velocity 



518.144 ftV volt “ aec 
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FIGURE H3 

Station 3 X-Array MS Fluctuating Axial Velocity 



153*69 ft /volt-sec 




co 



( OTXS 410 A) n - X^tooiga ITJTXV Sux^Bn^onx^ SW 



87 



FIGURE kk 

tation I 4 . X-Array MS Fluctuating Axial Velocity 



178.98 ft* /volt-soc 
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Station 6 X-Array MS Fluctuating Axial Velocity 
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FIGURE k7 

Station 7 X-Array MS Fluctuating Axial Velocity 




174. li -5 ft'/volt- 
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FIGURE 14-9 

Station. 9 X- Array MS Fluctuating Axial Velocity 
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FIGURE 50 

Station 10 X-Array MS Fluctuating Axial Velocity 



(3 = 525.87 ft /volt-sec 
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Station 2 X-Array MS Fluctuating Lateral Velocity 



315. 1+3 ftVvolt- 
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FIGURE 53 

Station 3 X-Array MS Fluctuating Lateral Velocity 



153*69 f tVvolt-sec 
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FIGURE 54 

Station 4 X-Array MS Fluctuating Lateral Velocity 
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FIGURE 55 

Station 5 X-Array MS Fluctuating Lateral Velocity 
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FIGURE 56 

Station 6 X-Array MS Fluctuating Lateral Velocity 



P Z = 161.12 ft '/volt- 
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FIGURE 57 

Station 7 X-Array MS Fluctuating Lateral Velocity 
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figure; 58 

Station 8 X-Array MS Fluctuating Lateral Velocity 




165.29 ft /volt-sec 
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Traverse Position (inches) 

FIGURE 59 

Station 9 X-Array MS Fluctuating Lateral Velocity 
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Station 10 X-Array MS Fluctuating Lateral Velocity 
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Station 1 X-Array MS Fluctuating Velocity Sura 



518.44 £t 2 /volt - 
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Station 2 X-Array KS Fluctuating Velocity Sum 
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(3^= 153-69 ft Vvolt-sec^ 
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Traverse Position (inches) 

FIGURE 6k 

Station 4 X-Array MS Fluctuating Velocity Sum 





178.93 ftVvolt-sec 
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Station 5 X-Array XS Fluctuating Velocity Sura 
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Station 6 X-Array MS Fluctuating Velocity Sum 
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Station 7 X-Array MS Fluctuating Velocity Sura 



171;. 1|5 ft ^volt-sec 
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165.29 ftV v °lt-sec 
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FIGURE 69 

Station 9 X-Array MS Fluctuating Velocity Sum 
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FIGURE 70 

Station 10 X-Array MS Fluctuating Velocity Sum 
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FIGURE 71 

Station 1 X-Array MS Fluctuating Velocity Difference 



P = 518*44 f t 2 /volt-sec' 
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p z = 315. Ij. 3 ft a /volt- 
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FIGURE 73 

Station 3 X-Array KS Fluctuating Velocity Difference 



53*69 ft '/volt-sec 
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FIGURE 7k 

Station I}. X-Array MS Fluctuating Velocity Difference 
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FIGURE 75 

Station 5 X-Array MS Fluctuating Velocity Difference 
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FIGURE 76 

Station 6 X-Array MS Fluctuating Velocity Difference 
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FIGURE 77 

Station 7 X-Array MS Fluctuating Velocity Difference 
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165.29 f t^/volt-sec 
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FIGURE 79 

Station 9 X-Array MS Fluctuating Velocity Difference 
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FIGURE 



Si 



The wind Tunnel Inlet 
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FidU HE 82 

The Wind l'unnel Test Section 
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PJdUHE 83 



Tho wind Tunnel Blower 
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Pi (Hi R i 8Ii 

i’he Test Section Traverse 
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FIGURE 85 

The Test Section Access Panels and 
Static Pressure Ports 
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F 1 uU Pul 66 



The Manometer Board 



129 




FIGUrik 87 



The Wake Generator 
With Abrasive Strip 
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FIGURE 88 



The Data Sensors 
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FIGURE 89 

The Test Section During Phase I 



132 






FIGURE 90 



X- Array and Normal Wire 
Anemometers 
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FIGURE 91 



The Anemometer Modules 
and Correlator 
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FIGURE 93 



Station 2 Turbulent Stress Comparison 
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FIGURE 91+ 

Station 3 Turbulent Stress Coraparison 
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FIGURE 95 

Station ij. Turbulent Stress Comparison 
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FIGURE 96 

Station 5 Turbulent Stress Comparison 
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FIGURE 97 

Station 6 Turbulent Stress Comparison 



APPENDIX A 



DEVELOPMENT OF THE GOVERNING EQUATIONS 
OF MOTION FOR TURBULENT INCOMPRESSIBLE FLOW 

First, the following equations are the basic Navier- 
Stokes equation for incompressible Newtonian flow: 
continuity : 



V • V = 0 



(A-l) 



momentum: 



U + V • V V = - - VP + v V 2 V 

3t p 



( A-2 ) 



Expanding the DEL terms they are 



3u , 3v , 3w _ A 

37 + 37 + 37 ' 0 



(A- 3) 



3U 4. ,, 3U 4. 4. ,, 9U 

ttt + Ur- + Vr— + Wr- = 

3t 3x 3y 3z 



1 3P L ,3 2 u j 3 2 u l 3 2 u> t A (l >, 

"n 77 + v( 2 + — 2 — 2 ' (A_I) 

P 9X 3 x 2 3 y 2 3 z 2 



3v , 3v , 3v , 3v 

rr + Ur- + Vr— + Wrr 

ofc 9 x 9 y 9 z 



1 3P , ,3 2 v , 3 2 v , 3 2 v 

p 3y ~ 2 .__2 ~ _2 



) (A-5) 



3x 



3y 3z 



3w , 3w , 3w , 3w 

rr + Ur— + Vr— + Wr- 

9t 9x 9y 9 z 



13P , ^ 3 w j 3 w , 3W\ t A 
-- r— + V r- + ~ + ~) (A-6) 

p 92 3x 2 3y 2 3z 2 



1*11 



Now substitute equation (A-l) into (A- i l) J (A-5)> and ( A— 6 ) 
to get 



au 
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3 2 
3x u 
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uv 

3y 
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3P 
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, 3 2 u 

v(— 2 

3x 
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3 2 u 
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CVJ 

CO 


3 1 
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3x 
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3y 2 




3z 2 


3v 
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|uv 

3x 
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3 2 

3y v 


X 




_l 


3P 
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/ 3 2 v 
v ( ? 

3x^ 
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3 2 u 
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3 2 v 


3t 
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~p 


3y 
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3y 2 


T 


3 z 2 


3w 
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3 

3x uw 


X 


1" w 
3y 
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_i 


3P 


X 


/3 2 w 
v( ? 

3x^ 


X 


3 2 w 


X 


3 2 w 


3t 


l 


i 


1 


p 


3 z 


1 


T 


3y 2 


1 


3 z 2 



With p and v constant and using the following definitions 

u = u + u 1 w = w + w ' 

v = v + v 1 P = P + P 1 

equation ( A— 3 ) becomes 

| x (u + u') + |- (v + v’) + |- z (w + w') = 0 ( A— 7 ) 



Since the integral over time of any fluctuating quantity is 
zeroj taking the time average of equation (A-7) results in 



3u 3v 3w _ 

3x 3y 3z 



(A-8) 



Then subtract this from (A-7) to get 



3u 1 , 3 v ' 3w ' _ 

3x 3y 3 z 



(A-9) 



1^2 



Now substitute the definitions above into equations (A-4), 
then 



~ t (u+u') + |- x (u+u ') 2 + (u+u' ) (v+v' ) + -| z (u+u' ) (w+w' ) 



cL <L d_ 

- " + v ^~2^ u+u '^ + ~2^ +u '^ + ~2^ +u '^ 

p ax ay az 



( A-10 ) 



and after taking the time average 



0 U. , 0 . 0 .<2.0 . 0 7 7” 0 . 0 \ 

5t + 9x u + 3x u + 3y uv + 3y u v + 3z uw + 3z u " <A - n) 



_ 1 IP , r 3 u j. 9 2 u ^ 3 2 «t 

- -p n + v[ 7 V 7 T r - ?- 1 

y ax ay az 



In a similar manner, one can get 



— + f-uv + u'v' + — v 2 + — v 1 2 + f-vw + v'w’ (A- 12 ) 

at ax ax ay ay az az 



2 — 2 — 2 — 

i ap L r a v L a v , a v-, 

= -“ X77 + VL o + o + oJ 

p 9y ax 2 ay 2 az 2 



7777 + I- UW + r- u'w' + I-..VW + 



at ax 



ax 



. .... 7 - V’W' + 7 T W 2 + 7- W ' 2 (A-13) 
ay ay az az 



— 2 2 2 

i ap L r- a w a w x a w-i 

- -7 V 7 + VL J + o + ?J 

p 9 ax 2 ay 2 az 2 



In the special case for a two-dimensional mean flow with a 
steady mean, equations (A-ll) and (A-12) can be written 
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■ip— + 1- uv 
ox 9y 



19,- ,2s 

T9x (p+pU } 



~ ~2- .2- 

L(^) + v[^-| + U] 

3y 9x 2 9y 2 



( A— 1 If ) 



^ v 2 
9x 9y 



■Hy (?+pv ' 2) 



^ ^2— r,2 — 

!x™ + v Cp + f-£ 

9x 9y 



(A-15) 



Then applying the standard boundary layer assumptions 



P = P ( x ) 



, p du 

from -p = -p U 
dx K 



e dx 



(A-16) 



and in addition make the following approximations 



P u 



• 2 « 



(A-17) 



« iiTHF 

9x 2 9y 2 3y 



(A-18) 



When these assumptions are substituted into (A-l4), the 
following equation results: 





1 9P. 
p 9x 



9_ 

3y 



U* V 1 



(A-19) 



Which can be manipulated into the form in equation (A-20) 
using continuity 



9u 
9 




1 9P 
p 9x 



9_ 

9y 



U * V ’ 



(A-20) 
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applying the Boussinesq laminarizat ion , the turbulent shear 
becomes 



9_ 

sy 



u' v' 



hy ' F & 



(A-21) 



then, T = -p u'v’, is the Reynolds Stress, and y^(x), the 
Eddy Viscosity, is a function of the flow field. If (A-l6) 
and (A-21) are substituted into (A-20) the result is 



u 



9u 

9x 



+ 




U 



dU 

( 

e dx 



v T (x) 



9 2 u 

ay 2 



( A-22 ) 



An alternate approach to solving these equations is 
possible without the boundary layer assumptions (A-l6, 17, 

& 18). Since there are three equations, (A- 8) , ( A— 1^1 ) , and 
(A-15) and six unknowns u, v, p, u* , v' , and u'v', three 
additional equations are required. By subtracting the mean 
equations from the full equations as shown in (A-10), three 
additional momentum equations for the fluctuating portion 
of the flov; result; however, several additional unknowns 
also appear: 



u t 3 , v i 3 ^ an( j u ' v * ^ . 

One might also formulate additional equations for these 
variables, but higher order correlation unknowns would 
appear in them, and additional equations would be required 



ll\5 



ad infinitum. It is customary at this point then to approxi- 



mate 

and 



2 2 

v' , u' , and u'v' by some function of the mean flow 
solve the system of equations. 
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APPENDIX B 



DEVELOPMENT OF THE EQUATIONS FOR 



THE HOT-WIRE ANEMOMETER 



I. Single Normal Hot-Wire 

First assume the following definitions 
e = e" + e ' 
u = u + u ' 

Then, since the linearized output of the hot-wire is 
proportional to u 



Setting averages equal to averages, since they must be 
independent of the fluctuations, 



e = C u 



(B-l ) 



or 



e+e 1 = C (u + u ' ) 



(B-2 ) 



e = C u 



( B— 3 ) 



or solving for u 




e 



e 



e 



(B-n ) 
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and 




Now setting the fluctuations equal 



(B-5) 



e ' = C u ’ 

u ' “ h e ' 




which from (B-4) becomes 




(B-6 ) 



( B— 7 ) 



(B-8) 



II. X-array Hot-Wire 

Now assume in addition, that 



v = v + v* 



e A = e A + e A ' e B = e B + e B ' 



(B-9 ) 
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then since the wires are separated by 90° 




ucos0 + vsinO 



Now note that 0 = 45° so 

e A = C A ( U + v) 
eB = C B ( u - v) 



(B-10 ) 



At the centerline of the flow the lateral velocity is 
approximately zero, so e^ = e R and, therefore, C A = Cg = C. 
For this reason, once the probe is adjusted to force e^ = e^ 
at the centerline the following relations hold 



e A = C (u + v) 
e B = C (u - v) 



(B-ll) 
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But if these equations are summed and solved for u 



u = 



_ 1 _ 

2C 



[e 



A 



e B^ 



(B-12 ) 



and similarly 

v = 2 c ^- e A " e B^ (B-13) 

Then if 3 s 1/2 C, and the definitions of the parameters are 
substituted into (B-12) & (B-13), the resulting equations 
are 



u + u’ = g[e A + e A ' + e B + e B ’ ] 



y + V - 6[e A + e A ’ - e B - e B ’] 



As in B-I, set averages equal to averages and 



u = $[e A + e B ] = 3 [e A + eg] 



v = 6[e A - e B ] = S[e A - e B ] 



setting fluctuations equal to fluctuations 



u’ = 3[e A * + e B ’] = 3 C ( e A + eg)'] 



v’ = 3[e A ’ - e B ’] = 3 [ ( e A - eg)'*-] 



(B-14) 

(B-15) 



(B-16) 



(B-17) 



(B-18) 

(B-19) 
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Then square and average to get 



u’ 2 = 3 2 [(e A + e B )' 2 ] (B-20) 

v' 2 = B 2 C (e A - e B )’ 2 ] (B-21) 

Returning to equations (B-ll) and setting primes equal to 
primes 



e' - C(u’ + v' ) 
e B ’ = C (u ' - v 1 ) 

square and average 

e A 1 2 = C 2 (u' 2 + 2u ' v 1 + v' 2 ) 
e B ' 2 = C 2 (u' 2 - 21?^ + v' 2 ) 

then subtract 

e.’ 2 - e ' 2 = ^uW 7 " 

A .d 



or 



/ 



U* V* 





] 



(B-22 ) 



(B-2 3 ) 



(B-24) 



(B-25) 



Now add equations (B-23) 



and 



e A ’ 2 + e B ' 2 = 2C 2 [u' 2 + v' 2 ] 



or 



u ' 2 + V ’ 2 = 26 2 [o a ' 2 + e B ' 2 ] 



(B-23) 



Then sum (B-20) and (B-21) 



u' 2 + v’ 2 = 6 2 [(e A + e B )' 2 + (e A - e B )' 2 ] (B-2^4 ) 



equating (B-23) and (B-2^), and rearranging 



so 



2[e A ’ 2 + e B ' 2 ] 



(e A + e B )' 2 + <e A - Og)' 2 



= 1 



(B-25) 



Where 



or 



P Te - P e ‘ l/ 2 P[u e 2 + v e 2 ] 



Pm ~ P 



2 C ^ ^ 2C ^ ^ e A + e B^e + ^ 2C ^ ^ e A e B^e 

P 



e 2 = 2( P Te ~ P e ) 



p [(e A + e B^ + ^ e A e B^ -'e 



1 r~ 



and 



2 ( P Te P e ) 
P 



is a known array of data from the pitot-static phase of the 
study. 
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APPENDIX C 



DEVELOPMENT OF AN INTEGRATED FORM OF THE 
MOMENTUM EQUATION FOR THE CALCULATION OF 
THE LOCAL TURBULENT SHEAR COEFFICIENT 

First, the basic equations of motion for a two-dimensional 
mean flow are 

continuity : 



x y 



(C-l) 



momentum: 



uu + vu 
x y 



— P + — T 

P x p y 



( C-2 ) 



then integrate equation (C-l) from zero to a point y 



u 



v 



X 



y 




y 

dy = -/ v dy 
0 y 



y 

/ u dy = - [v + f ( x ) 3 
0 x 



y 



0 



thus 



y 

v(y) = - / u x dy + v 



0 



15 ^ 



is zero, so 



But 



v 



y=o 



y 

v(y ) = - / ' u dy (C-3) 

0 x 



Now simply and multiply equation (C-l) by u and add to 
(C-2) to get 



2uu + (vu + uv ) = - — P + — (C-4) 
x y y p x p 

where the terms on the left can be recognized as differentials 
of the form 



2uu x = (u 2 )x 

VU y + UV y = ( UV >y 

thus equation (C-^) becomes 

j/ 

P 

and, since v(0) = 0, when this equation is integrated from 
y = 0 to an arbitrary point y = y , 

y y y 

/ 2uu dy + / T (uv) dy = / T 2uu dy + u(y )v(y ) 
0 x 0 y 0 tt 



2uu x + (uv) y 



- - P + 
P x 



and the result, when equation (C-3) is substituted for 
v(y ), is the following relation 
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y T y T y TP 

/ 2uu dy - u(y ) / u dy = - / ~dy 
0 x T 0 x 0 p 



T(y ) 

+ — 

P 



Which can be rearranged into the following form, 

y x y x 

/ u [2u-u(y )]dy = - i / P dy + i T (y ) (C-5) 

0 x T PqXpt 

And can be solved for the coefficient of shear at y 

x 

(C T = T) 

pu e 

y T y T 

C(y T ) = — V" f p x dy + T" f u [2u-u(y )]dy] 
T T pu e 0 u/ 0 x T 



An alternate approach to this development would be to 
substitute equation ( C— 3 ) into (C-2) 



— 3u 3u r y / 9u \ _ 1 3P , 1 3 t 

U 37 ’ 3y </ ( 37 )dy " " p 8x p Sy 



(c— 7 ) 



Then non-dimensionalize by division by the freestream 
velocity squared, and use the product rule of the calculus 
to get the following form of the momentum equation: 



+ (-^) 2 - (-H-) 7 (-±) dy-(-a-) <^) 7 (— — )dy 

u e u e x u e Fe u e y 0 u e x u e y u e 0 u e 



-S- P x + (-f > 

p u 2 



(C-8 ) 
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2 T 

V/ here C = p > and the subscripts x, and y indicate 

P U e 

partial differentiation with respect to x, and y respectively. 
Now, if the condition that equation ( C— 8 ) be satisfied at 
every point in the flow is relaxed to allow it to be satis- 
fied on the average between y = 0 and an arbitrary point 
y , then ( C— 8 ) can be integrated by parts to get the 
following relation for 



W 



y x u(y ) u 

= / [2(^-) —IK—) +(—)( — 

0 u e u e u e x u e u e 



)]dy 



+ 



0 



y 

/ 



T 1 

2pu e 2 



P 



x 



dy 



( c— 9 ) 



Now if the original definition of the shear term 



t = - p u ' v* 

is substituted back into equations ( C— 5 ) and (C-9), the 
following relations result: 

1 

y T y T 

u’v' = / ulu(y ) - 2u]dy - ^ ^ P dy (C-10) 

0 x t p Q x 

2 (*j— ) ] C +(-j L )(^)]dy 

e e x e e 

(C-ll) 



P t u (y ) 

U * V ' = u / f [ 

0 



xP x 

V ^ dy 
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Finally , to check the validity of these equations, let 
y^ go to infinity in the limit and assume that 



P = - p u u 
x e e 



x 



Then, equation (C-5) becomes 



2 / uu dy - u / u dy - / u u dy = 0 

0 x e p , 0 e x 



Since the shear, and therefore C , is zero at °° and then 

3 t 5 

OO 2 00 OO 

2 / (^-) dy - / [(u u) - uu ]dy - / u u dy=0 

0 ^ x 0 x e x 0 ee x 



or 



/ C(u ) - (u u) + uu - u u ] dy = 0 

0 x e x e x e x 



thus, combining terms 



, r 2 U / U 

' [u e V ( — 
0 e e 



1,)] dy + / u u (^- - 1 )dy = 0 

x 0 e e x u e 

( C-12 ) 



and if the flowing parameters are defined 



6 V t ( t- 1)dy; 



6 * = / (^- - 1) dy 



u 



0 u e 
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then equation (C-12) may be written 



[u 8] + 



u u 
e e. 



0 



which when the first term is differentiated becomes 



2 d0 , „ a , _ 

u -r— + 2uu 0 + uu o — 0 

e dx e e e e 

X X 



or dividing by u 



g + _!i [2e + 6*] - 0 



(C-13) 



which is simply the well-known momentum integral equation 
for a boundary layer. A similar procedure will reduce 
equation (C-9) to (C-13) as well. 
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APPENDIX D 



DATA CHECKS 

In an experimental study, it is very important to have 
some idea of the quality of the data that results from the 
experiment. In this Appendix, several methods of checking 
the accuracy of the data in this study are discussed. 

CONTINUITY CHECKS 

Using the continuity equation, it is apparent that, 
since the mass flow must be constant, 
wall 

/ udy = constant (D-l) 

0 

Therefore, this equation was added to the program described 
in Appendix E. The results of this check for the pitot- 
static data varied as much as 30 % over the entire test 
section, yet, if the largest and smallest value are neglected 
the range is narrowed to ± 10 %. The author believes that 
much of this difference is due to the inability to position 
the pitot-static probe near the wall in some cases, while 
in others, the head of the probe could be placed adjacent 
to the wall surface. This problem is especially critical 
at Station 1 where the boundary layer is thin, and mass 
flow near the wall, at free stream velocity, comprises a 
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significant part of the total mass flow in the tunnel. 

See the computer output in Appendix P. 

MOMENTUM CHECKS 

Recall the integral momentum equation derived in 
Appendix C, 

U e 

If + — - [20 + 6*] = 0 (C-13) 

u 

e 

which is a special case for the. free stream where C is 
zero. Taking this equation as 



El 



d_0 

dx 



E2 = 



du 



u dx 
e 



E3 = 20 + 6 



then the non-dimensional residual of this equation was 
programmed as 



Epsilon = (El + E2 * E3)/(E1 - E2 * E3) (D-2) 



and its values are listed in the computer output. 

In general, this ratio showed the data to correlate 
fairly well with the restrictions imposed by the momentum 
integral equation with an average residual of less than 



l6l 



25 % and as little as + 15 % for most stations; however, 
Station 7 data was completely out of line with a residual 
greater than two. The primary reason for the failure of 
Station 7 to check satisfactorily was an inclusion of a 
significant portion of data at the upper portion of the 
traverse that reflected boundary layer interference. The 
inclusion of the resulting velocity decay in the calcula- 
tion of the momentum, and displacement thicknesses created 
a large error in their values. This error is almost 
certainly the cause of slightly high values of Station 8 
as well. 



Having calculated the local turbulent shear coefficient, 
it was possible to also integrate to the residual in the 
energy integral equation. 



Therefore taking the term on the left as A, and that on 
the right as B, the relation was programmed and. 



ENERGY CHECKS 




(D-3) 



where the energy thickness <5^ is defined as 



6 




( D— 4 ) 



ZETA = (A - B)/( A + B) 



(D-4) 
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The results of this data check were inconclusive. The 
magnitude of one of two terms above (A and B) was so much 
larger than the other that the difference over the sum was 
due solely to an error in the programming vice an inaccuracy 
in the raw data. 



HOT-WIRE CHECKS 

In Appendix B it was shown in equation (B-25) that these 
terms were set equal to n, and values for n are listed in 
Table I. It was immediately obvious in this case that the 
error, although not large, was directly related to the 
strength of the turbulence in the flow. Thus, the magni- 
tude of the lack of algebraic correlation in this check 
may have been due to an inability of the equipment to 
handle the rapid variations caused by severe turbulence. 

It is important to note here, that a change of as little 
as 10$ in the individual values of the parameters used 
to calculate n would bring even the worst value of n back 
to within 10$ of 1.0. The author felt that the data from 
Phase III of the experiment was well within the bounds of 
acceptable experimental data spread. 
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APPENDIX E 



DATA REDUCTION COMPUTER PROGRAM 



A. GENERAL 

The Fortran IV computer program included in this Appen- 
dix is a revised form of that developed by the study of 
Reference 1 . The program inputs the pitot-static data of 
phase I of the experiment, and a limited amount of data 
from phase II to calculate the wake momentum thickness, 
displacement thickness, energy thickness, local turbulent 
shear coefficient, and Eddy Viscosity. 

At the inception of the present study the portion of 
the program dealing with the turbulent shear coefficient 
did not produce accurate values for C . The author of 
Reference 1 was unable to find a reason for this error 
since the equations appeared correct, and there were no 
obvious coding errors. He postulated that the unreasonable 
values were a result of an inability of a forward finite 
difference scheme to accurately approximate the partial 
derivatives in equation ( C— 9 ) of Appendix C. Since this 
appeared to be a reasonable theory, this author attempted 
to obtain values for the derivatives by two additional 
methods : 

1. A least squares polynomial approximation was fitted 
to the data at a given position y from the 
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centerline then this polynomial was differentiated. 

2. A central difference scheme. 

Unfortunately, these improvements resulted in no signi- 
ficant change in the values of x. Thus, the values for 
the various derivatives were printed out in detail, and 
an order of magnitude analysis was carried out. At this 
point, it was not difficult to see that the values of the 
parameter U were exceptionally smooth in the x-direction, 
and obviously not a problem to approximate, yet the values 
calculated by the IBM-360 were not even similar to those 
calculated by hand. The author then attempted to follow 
the instructions programmed for the computer and quickly 
discovered that a subtle error was hidden in the program. 
Since each station had a different number of points, the 
index corresponding to a point y above the centerline was 
different at each station; therefore the arrays requiring 
differentiation were read into a dummy variable where 
J = 1 corresponded to the centerline for each station. 

This correction produced values for the local shear coeffi- 
cient that were of the right order of magnitude, and the 

I 

central difference scheme was considered to produce the 
best values for the derivative. Some of the hot-wire data 
was also analyzed using this program, but the values for 

were not very similar to those calculated from the pitot- 
static data. Since the pressure term was found to be the 
driving force of the equation for C during the order of 
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magnitude analysis, the author felt that approximating 
the static pressure data with the data from phase one would 
be unlikely to produce realistic values for the parameters, 
and no additional attempt was made to use the IBM-360 to 
reduce the hot-wire data. 

B. NOMENCLATURE 

The following parameters and arrays were used in the 
program; their respective dimensions are shown in 
parenthesis : 



ALV 


array used to store interpolated values of 
dynamic pressures (feet) 


AMI XL 


local Prandtl mixing length in the wake (fee) 


CTAU 


local turbulent shear coefficient for each point 
in the wake 


DELRAT 


ratio of succeeding DELSTR' s to that at the first 
data station 


DELSTR 


wake displacement thickness (feet) 


DELX 


axial distance between succeeding probe tip 
positions (feet) 


DIAMTR 


diameter of wake generator (inches) 


EDDY 


eddy viscosity 


EPSIL 


residual in momentum integral equation when 
experimental values for THETA and SHAPE are 
substituted into it 


LABEL 


label for subroutine DRAW 


MN 


number of data points laterally through the wake 


MX 


number of axial data positions 


MY 


number of points (always odd) into which the 
wake is divided 
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p 


array of local static pressures for each point 
in the wake (psf) 


PATM 


atmospheric pressure (inches of mercury) 


Q 


freestream dynamic pressure (psf) 


SHAPE 


wake shape factor (DELSTR/THETA ) 


THETA 


wake momentum thickness (feet) 


THETAP 


inclination of static pressure manometer 
(degrees ) 


THETAV 


inclination of dynamic pressure manometer 
(degrees ) 


THRAT 


ratio of succeeding THETA's to that at the first 
data station 


TITLE1 


title for subroutine DRAW 


TITLED 


title for subroutine DRAW 


U 


array of velocities for each point in the wake 
(feet per second) 


URATIO 


ratio of local wake velocity to local freestream 
velocity 


XIN 


array of probe traverse positions through the 
wake (inches) 


XOUT 


interpolated probe traverse positions for each 
of the MY points in the wake (inches) 


XPOS 


probe tip position measured from downstream 
side of the wake generator (inches) 


YOUT 


interpolated manometer readings (for YQ and YPS) 
for each of the MY points through the wake 
(inches ) 


YPS 


static pressure manometer slant distance height 
difference (inches) 


YQ 


dynamic pressure manometer slant distance 
height difference (inches) 


MU( X) 


kinematic viscosity 


MM 


index of the maximum velocity 
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MT 


index of the maximum velocity in the upper 
half of the data 


T 


dummy array of U for differential 


TRATIO 


dummy array of URATIO for differential 


ADP 


dummy array of ALP for differential 


D 


dummy array of for calculations Eddy Viscosity 


NUM 


array of MN 


MBM 


array of MM 


MAD 


array of MID 


MUT 


array of MT 


MID 


index of centerline 



C. SUBROUTINES 



INTEG 


uses Simpson's rule to numerically integrate the 
function F. It must be supplied with K points 
(always odd) equally spaced, DELX units apart 


QTFE 


SSP3 library subroutine for integration using 
trapazoidal rule 


QSF 


SSP3 library subroutine for integration using 
Simpson's Rule 


FILL 


used to fill in arrays to maximum dimension with 
freestream values 


PREIN 


determines midpoint of asignally dimensioned 
array 


UTOP 


determines MT, & MM 


SPLIT 


takes the entire array and puts the values above 
the centerline into a second array in the 
input parameters 


D. USAGE 





The details on the use of this program are listed in the 
comments at the head of the program. The deck requires 120K 
to compile and 230K to run. 
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APPENDIX E 



DATA REDUCTION COMPUTER PROGRAM 

^ ^ ^ ^ ^ A ^ ^ ^ A ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ * 4? ^ sic * 

A HOT-WIRE ANEMOMETER STUDY OF FREE TURBULENT MIXING 
IN AXIAL PRESSURE GRADIENTS. 

#{- fc sjc $c 7 %. % & Jfc ### # 'Jr.7z%. sgs$$#al : & 

THIS PROGRAM WAS DESIGNED TO REDUCE P I TOT- STA T I C , 
AND HOT-WIRE ANEMOMETER DATA TO DETAILED DATA 
CONCERNING THE PHYSICAL PARAMETERS OF THE AIR FLOW. 

IT WAS SPECIFICALLY DESIGNED FOR USE IN THE THESIS 
PROJECT OF LT. M. E. DEARNEY IN 1972, HOWEVER, THE 
PROGRAM WILL BE USEFUL TO ANYONE WHO HAS DATA OF A 
SIMILAR NATURE. 

TYPE OF DATA REQUIRED 

EITHER PITOT STATIC, OR HOT WIRE DATA (WITH A 
VELOCITY REF.) AT EQUAL INTERVALS ACROSS A PARALLEL 
FLOW 2— DI MENS I ONAL WIND TUNNEL. THE PROGRAM IS 
PRESENTLY SET UP TO ANALYZE 3 SETS OF DATA OF 12 
STATIONS EACH, THE FIRST BEING PITOT STATIC VALUES 
AND THE SECOND AND THIRD, HOT-WIRE DATA. ONLY MINOR 
MODIFICATIONS WOULD BE REQUIRED TO EXTEND THIS TO 
MANY MORE EXPERIMENTS. 

I NPUT: 

1. SIX TITLE DESCRIPTION CARDS 

2. ONE LABEL, 1 NTERGERS FROM 1-12 

3. FIRST RUN TITLE 

-4. FIRST 2 LINES, NUMBER OF STATIONS 
DIAMETER OF C YCLI NDERl NOT REQUIRED) 

5. ONE OR MORE CARDS WITH MANOMETER ANGLES 

6. ONE OR MORE CARDS WITH ATMOSPHERIC 
PRESSURES. 

7. TITLE OF FIRST STATION, FIRST RUN 

8. FIRST 3 LINES, INTERGER NUMBER OF POINTS 
I <17 5 . 

9. REFERENCE VALUE OF FIRST POINT 

10. DYNAMIC PPESSURES IN I NC Hr S ( PT- PS ) 

11. STATIC PRESSURES IN INCHES(PS) 

12. SAME DATA 7-11 FOR FOLLOWING STATIONS 
11. SECOND RUN TITLE 

14. NUMBER OF STATIONS 
DIAMETER OF CYCLINDER 

15. SAME AS 5,6, L 1 

16. NUMBER OF POINTS 

17. REFERENCE VALUE 

18 . VOLTAGE DATA 

19. THEN CONTINUE FOR THE REST OF THE STATIC 
THE THIRD RUN IS THE SAME AS THE FIRST, 

MAJOR VARIABLES 

U-VELOCITY 

CTAU-COEFFI Cl ENT OF STRESS 
MU( X) , VI SCOUSI TY MODEL 

p__p D F Q C ! I P F 

UR A T 10— U/U (FREE STREAM) 

AL 3 - P ST AT I C / 12 

ALV- PDYNAMIC/12 

EDDY- EDDY VISCOUSITY 

AMI SL- REYNOLD'S MIXING LENGTH 

MAJOR SUBROUTINES 

INTEG- INTEGRATE BY SIMPSONS RULE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

cc 

c 

c 

c 

c 

c 

c 

c 

c 

50 

51 

52 

53 

54 

55 

56 

57 

58 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 



QTFE- INTEGRATE BY TRAPAZOI DAL RULE 
FILL- FILL THE VELOCITY AMD PRESSURE ARRAYS 
FOR COMPUTATIONAL PURPOSES* 

PREIN- DETERMINE THE MIDPOINT OF THE WAKE 
DEFECT . 

SM- USED IN CONJUNCTION WITH LEAST SQUARES 
ROUTINE FOR TAKING DERIVITIVES, 

1-U, 2= DU. 

SPLIT- DETERMINE THE MIDPOINT OF THE FLOW 
UEE- DETERMINE THE MAX FLOW VELOCITY 






IMPLICIT REAL*4 (A-H,0-Z) 
REALMS TITLEK12) 

REALMS TITLE2 (12) 

REAL*8 T I TL E4 ( 12) 

REALMS L A 3 EL 1 / ' 

REAL*8 LABEL?/' 

RE AL*8 LABEL3/' 

RE AL--4 LABEL (12) 



DELS '/ 
THET' / 
SHAPV 



DIMENSION YFAT! 12) , YFI T( 12) , YFOT! 12) , YFUT! 12 ) 

DIMENSION ALP( 12.175) * ALV{ 12 , 175) »U( 12 , 1 7 5 ) , P ( 12, 175) 
DIMENSION AMU( 12 ) , MUT( 12) 

DIMENSION A(12,175),B(12,175),C(12,175),D!12,175) 
DIMENSION EDO/ (12,175) , AM I XL (12, 175),X0UT( 12,175) 
DIMENSION URATI 0(12, 175) , C TAU ( 1 2 , 1 75 ) , YQ ( 1 75 I ,YPS(175) 
DIMENSION UR( 1 75) , UP ( 175) , UC( 17 5) ,UD(17 5) ,F( 17 5) 
DIMENSION YOUT ( 175 ) ,XI N( 175 ) 

DIMENSION A X I S ( 175) 

DIMENSION DELS T R ( 1 2 ) , THETA ( 12) , SHAPE C 12) , EPSIL ( 12 ) 
DIMENSION T(12 ,175) ,ADP(12, 175) 

DIMENSION 0(12) , DE L RAT ( 12) , THR AT ( 1 2 ) , THETA V( 12) 

DI ME NS I ON THETAP (1 2 ) , PATM( 12 ) , XPOS ( 12 ) , DEL Y { 12 ) 

DIME NS I ON WX( 12) , WH ( 12) , I W( 12) 

DIMENSION XF IT ( 12) 

Dl ME NS I ON RAKE (175 ) 

DIMENSION FAKE ( 175) t MAD { 12) ,MBM( 12) 

DIMFNS ION NUM( 12 ) 

COMMON /A/ H 

DATA GAMA/62. 32/, RHO/O. 002378/ 

H = 0 . 1 E 0 
Y=H/12 . 0 



SET INPUT/ OUT 5 UT FORMATS 

FORMAT ( ' 

1 ' ) 

FORMAT ( 12, F10. 2 ) 

FORMAT! 8F10. 1) 

FORMAT { 8F10. 2) 

FORMAT ( ' 

1 ' ) 

FORMAT ( 13) 

FORMAT ( 8F10.3) 

FORMAT (///////' ' , T20 , 6A 8) 

FORMAT (/////// • ' , T 20 , 12A4 ) 

FORMAT ( //////////' ' ,T33 ,' ECHO-CHECK DATA'///) 

FORMAT!//' ' , T 20 , ' THE TA V’ , ( ' • , T33 , 8F 1 0. 2 / /) ) 

FORMAT ( // ' ' ,T20, ' THETAP ',( ' ' , T33, 8F 1 0. 2/ / ) ) 

FORMAT!//' ' ,T20 , ' PATM' « ( ' ' ,T3 3 , 8 F 1 0 . 2 / / ) ) 

FORMAT!//' ' , T 20 , ' XPOS ' , ( ' ' , T33 , 8F 1 0. 2 / / ) ) 

FORMAT!///' ' ,T53, ' MN=« , 13///) 

FORMAT!//' ' ,T 10, ' XIN' , ( ' ' , T2 1 , 8 FI 0 . 3 / / ) ) 

FORMAT!//' ' , T 10 , ' YQ ',(' ' , T2 1 , 8F 1 0 . 3/ / ) ) 

FORMAT!//’ ' ,T10,' YPS' , (' ' , T2 1 , 8 F 1 0 . 3/ / ) ) 

FORMAT! ' 1' ,T50, '***RESULTS OF DATA REDUCTION***') 
FORMAT (' 1* ////// ' ',T50, 'AXIAL ST AT I ON ', I 2 ,' OUT OF', 12 

1 ' STATIONS' /////) 
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•nRMATl ' ' , T18, ’ D E L S T A R = * ,G12. 5,5X, * THE T A = 1 ,G12.5,5X, 

RM AT? // / / // 5 ' ’ 8 ,' 01 NT 7 ! T3 0 ] • U ' , 1 4 0 , ' U /U E • , T 7 0 , • P • 

^•0RMAT ( M 0^'t9,I3,T24,G12.5,T4A»G12.5,T64,G12.5,T84, 

: '=ORMAT( ' 1' .T't4» '***SUMMARY OF INTEGRATED WAKE PROPER', 

TT \ 

nRMAT /////* ‘ ,T3, ‘AXIAL* ,T75,' THETA*/* • , T2 , ' ST AT I ON 

• °T16,*X/D‘,T29, ' DELST AR ' , T 44, • DEL RAT 1 0 • , T 60, * THE TA * , 

‘-ORMATP 1' ?T49^'*^RPYN6LD"S 0 MixiNG 5 LENGTH DATA*** 1 ) 

FORMAT (//////• *,T 61, -STATION* ,13///' • ,T44,' POINT- , 

T59 , ' EDDY VISCOSITY* ,T79,' MIXING LENGTH*) 

"OR M ATI* * , T45 , I 3 * r60,G12.5,T80,G12«5) 

-ORMAT ( 6 A8 ) T , T „ , 

: ORMAT ( / / - * ,' L AS F = 1 ,13) 

: ORM AT ( 12A4 ) 

Format! * i* //// //«' a ^T8 L , C ' 1 point* 1?3o( 'u* ,T48 , * u/ue* ,T70, 
, T88 , * CTAU * ) 

\ -ORMAT ( ' 0* { T60 , ' TURBULENCE CHECKS* / / /T1 0 , * ST AT I ON • , 
15X, ' ET A ' , 12X, * MU( X ) * / / ) _ 

r FORMAT ( ’ 1 ,T14*I2,11X,F12.4,5X,F12.4/) 



READ(5 ,90 ) 
RE AD ( 5 , 90) 
RE AD ( 5 , 90 ) 
READ(5 ,92) 



T I TL El 
Tl TLE2 
Tl TLE4 
LABEL 



ZERO ALL ARRAYS BEFORE COMPUTATION BEGINS 



THIS LOOP CARRIES THE PROGRAM THRU THREE SETS OF DATA 






*#** ******** «*******$ ** 5 =* 



DO 1000 KK=1 , 3 
Z = 0. 0 

DO 1 1=1,12 
DELRAT ( I ) = Z 
DELSTR ( T ) =Z 
DEL Y ( I ) = Z 
EPS I L ( I )=Z 
PATM( I ) =Z 
0(I) = Z 
SHAPE ( I ) = 7. 

THETA ( I )=Z 
THET AP ( I ) = Z 
THETAVl I ) = Z 
THRAT ( I ) =Z 
AMU ( I ) = Z 
XPO$( I ) = Z 
DO 1 J = l, 175 
A( I » J ) = Z 
8 1 I , J) = Z 
C( I , J ) =z 
CTAU (I , J) = Z 
D ( I , J)=Z 
F ( J ) =Z 
P( I, J) = Z 
URAT I 0 ( I , J ) = Z 
XOUTl I , J) =Z 
GO TO (36, 1, 1) , KK 
U ALP ( I , J ) =Z 
ALV( I , J) =Z 
U(I, J)=Z 
CONT INUE 
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J 7=MZ 
MZ=MZ+1 

DO 864 J =MU v 1 7 5 
YPS(J)=ALP( I,ME)*12. 
YPS ( J7 ) =YP S ( MZ ) 
IF(J7.E0.1) DO TO 862 
J7 = J 7- 1 



864 CONTINUE 
862 CONTINUE 
GO TO 43 
41 RE AD ( 5 , 56 ) 



(YPS (J ) ,J=1 ,MN) 



ECHO-CHECK INPUT 



43 WRITE (6,95) 
WRITE! 6, 54) 
WRITE* 6 ,75 ) 
WRITE ( 6,76) 
WRITE{6,95) 
WRITE (6,77) 
WRITE! 6, 95) 
WRITE (6,78 ) 
M Y = M N 
NUM ( I )=MN 

CALCULATE 



MN 

( XI N( J ) ,J=1 , MN) 
( Y 0 ( J ) , J = 1 , MN ) 

( Y P S ( J ) , J = 1 , MN ) 



U VELOCITY 



IF(KK . GT »1 ) GO TO 3 7 
DO 3 J=1,MN 
YOUT! J )=YO( J ) 

ALV ( 1 , J ) =Y OUT ( J ) / I 2 , 0 

U( I , J 5 - SORT ( 2. 0*GAMA*SIN{THETAV{I ) )*ALV ( I , JI/RHO) 
) CONTINUE 
GO TO 38 
37 MM=MBM { I ) 

UDUM=U l I , MM ) 

CALL PREIN(YO, MN, MI D) 

MM=M I D 

DO 632 J=1,MY 
632 U ( I , J ) =YQ( J) 

UTOP = UE E < U, MN , UMAX, MM, I ,MT) 

MV =M BM ( 1 ) 

ADUM=AL V ( I , MV) 

DO 39 J = 1 , MY 

ALV ( I , J ) = ( A DU M *YQ ( J ) ) / YQ ( M M ) 

U( I , J) =( UDUM*YO( J) ) /YO(MM) 

39 CONTINUE 



DETERMINE MAXIMUM FREESTREAM VELOCITY 



38 CALL PREIN(YQ,MN,MID) 

MM = M I D 

UTO? = UE E ( U, MN , UMAX , MM , I , MT ) 

0( I ) = 0 » 5*RHO*U( I ,MM)**2 

WRITE(6,55) MID 

WRITE (6,55) MM 

WRITE16, 55) MT 

WRITE (6, 53) UMAX 

WRITE! 6,53) UTOP 

MAD ( I ) = M I D 

MB Ml I ) =MM 

MUT ( I ) =MT 

CONTINUITY CHECK 

DO 99 J = 1 , MY 
99 F ( J) =U( I ,J) 

J4= 2 

CALL SPLIT (F, FAKE, MT, MID, J4) 
CALL INTEG(FAKE,J4,Y, RESULT) 
WR ITE ( 6 , 93 ) RESULT 
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CALCULATE DISPLACEMENT THICKNESS 
DD 31 J= It MY 

URAT I G ( I * J) =U( I , J)/U( I ,MM) 

F ( J) =1. 0-URATI OH , J) 

31 CONTINUE 
J1 =2 

CALL SPLIT (F,FAKE»MY » M I D , J 1 ) 

CALL OTFE(Y»FAKE»RAKEtJl) 

DELSTR! 1 ) = RAK E ( J 1 ) 

DELRAT ( I ) =DELSTR( I ) /DELSTR(l) 

CALCULATE MOMENTUM THICKNESS 

DO 4 J = 1 » M Y 

F ( J ) =UR AT I 0 ( I t J )*( 1 .O-URATIO! I , J ) ) 

4 CONTINUE 
J 2 = 2 

CALL SPLIT (F, FAKE, MY, MID, J2) 

CALL QTFEiY, FAKE, RAKE, J2) 

THETA! I J=RAKE( J 2 ) 

THRAT! I ) =THETA ( I ) /THETA! 1 ) 

CALCULATE SHAPE FACTOR IN MOMENTUM INTEGRAL 

SHAPE! I )=DELSTR! I) /THETA(I ) 

CALCULATE THE PRESSURE FIELD 

DO 5 J = 1 , M Y 
I F ( KK . GT . 1 ) GO TO AO 
X=XOUT ( I ,J) 

YOUT! J)=YPS( J) 

ALP! I , J ) =Y OUT ( J ) / 1 2 . 0 
F ( J ) =UR A TI 0! I , J) 

40 P(I,J)=PATM( I)- GAMA*SIN( THETAP ( I ) )*ALP( I , J) 



5 CONTINUE 



CALL DRAW! MN,XIN,F , 0,0, LABEL ,TI TLE1 ,0,0 ,1 , 1,2 ,2,6 , 6,1 



6 CONTINUE 



FILL IN FUNCTION ARRAYS TO MAX DIMENSION FOR 
CALCULATIONS 



CALL FILL! NUM, 1 , U ) 

CALL FILL! NUM, 1 , ALP) , 

PREPARE DUMMY VARIABLES CENTERED ON THE FLOW 
CENTERLINE 

DO 711 LT B= 1 , 2 
IF(LTB.EQ.2) SO TO 712 
DO 431 1=1, MX 
MI D=MAD! I) 

DO 451 J= MID, 175 
JK= J-MI D+l 
T! I , JK ) =U ! I , J ) 

ADP! I , JK) =ALP( I , J) 

451 CONTINUE 
GO TO 714 
712 DO 713 1=1, MX 
MI D=MAD( I ) 

JL =MI D +1 
DO 713 J=M I D , 1 00 
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JK= J-M I D+l 
Jl - J L- 1 

T ( I , JK) = U( I , JL ) 

ADP ( I, JK)=ALP! I, JL ) 

IF { JL. E0.1 ) GO TO 7 IT 

713 CONTINUE 

714 CONTINUE 

CALCULATE RESIDUAL IN MOMENTUM INTEGRAL EQUATION 
I =1 

DEL X=XPOS( 1 + 1) — X D 0 S ( I ) 

MV=MBM ( 1+1 ) 

MM=M8M ( I ) 

E 1= ( THETA ( I + 1 ) - THE TA( I ) ) /DELX 
E2 = ( THET A ( 1 )*( SHAPE! I )+2. ) ) 

E 3 = ( U ( 1 + 1 » MV) - U ( I , MM)) /(DELX*U(1»MM) ) 

EPS IL ( I )=! E1+E2*E3)/!E1-E2*E3) 

DO 7 1=2, MX 1 

DELX = XPOS( I + 1 ) — XPOStI — 1) 

Ml =MBM ( 1-1 ) 

M2=MBM( I ) 

M3=MBM ( 1+1) 

MY = NUM ( I ) 

El = ( THETA! H-l) -THETA! 1-1 )) /DELX 
E2= THETA! I )*( SHAPE! I ) +2 . 0) 

E3 = (U( I +1 1 M3 ) -U ( I-1,M1 ) ) / ! DELX^U ( I, M2) ) 

EPSILd ) = ( E1 + E2 : ! < E3 ) / ( E 1 -E2* E3 ) 

7 CONTINUE 
I = MX 

DEL X=XPOS( I )-XPOS( 1-1 ) 

MV=MBM( I ) 

MM=MBM ( 1-1 ) 

E 1=( THETA! I ) -THETA! I— 1 ) )/DELX 
E2=( THETA! I )*( SHAPE (I )+?..) ) 

E3 = ( U ( I * MV ) — U ( I-1,MM))/(DELX*U( I, MV) ) 

EPSILtI )=!E1+E2*E3)/(E1-E2*E3) 

CALCULATE LOCAL TURBULENT SHEAR COEFFICIENT IN WAKE 
STATION ONE - FORWARD DIFFERENCE 



1 = 1 

MM=MBMl I ) 

MY=NUM( I ) — MAD ( I ) -t-1 
DELX = X PQS ( I + 1)-XP0S(I ) 

DO 576 K=1,MY 

DO 575 J = 1 , K 

A! I » J ) = 2.*T ( I » J )--T ( I , K ) 

B ( I < J ) = ( T ! I + 1 , J ) — T ( I , J ) ) / DELX 

C! I , J)=( ADP! 1+1, J)-ADP { I , J) ) /(DELX*ALV( I ,MM) ) 

575 CONTINUE 
J4 = 0 

DO 574 J = 1 , K 
J4= J4+1 

574 Ft J) = A C I ,J)*B! I , J) 

CALL QTFE! Y, F, RAKE , J4) 

J5 = 0 

DO 573 J=1 , K 
J 5= J 5+ 1 

573 F! J) =C< I , J ) 

CALL QTFE! Y,F, FAKE , J5) 

CTAU! I,K)=(2./U( I,HM)**2)*RAKE(J4)-FAKE( J5) 

576 CONTINUE 

STATION 2 TO ( MX— 1 ) CENTRAL DIFFERENCE 

DO 10 1=2, MX1 
DELX=XPOS ( 1+1) -XPOS! I ) 

MM=MBM! I ) 

MY=NUM ( I ) - MAD ( I ) +1 
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M I D= 1 

DO 10 K = MI D , MY 
DO 9 J=MID,K 
MA= J— MI D+l 

A! I , J) = 2.#T (I,J)-T(I,KJ 
B(I.J)=(T(I+1,J)-T(I-1,J))/(2.*DELX) 

C( I , J ) ={ ADP! 1 + 1 , J J-ADP ( 1-1 v J ) )/ ( 2.*DELX*ALV( I , MM) ) 

9 CONTINUE 
J4=0 

DO 701 J = M ID, < 

MA = J-Ml D+l 
J4= J4+ 1 

701 F(MA) = A( I, J)*B( I ,J ) 

CALL OTFE( V, F, RAKE , J4 ) 

J 5=0 

DO 702 J = M I D , < 

MA= J— MI D+l 
J 5 = J 5+ 1 

7 02 F ( M A ) = C ( I , J ) 

CALL QTFE(Y»F, FAKE, J5 ) 

CT AU ( I,K)=( 2./U( I,MK)**2)*RAKE{ J4)-FAKE( J5J 
10 CONT I N U E 

I F ( LTB . EO. 2 ) 30 TO 718 
DO 521 1=1, MX 
MY = NUM( I ) — MAD ( I ) +1 
MI D = MAD ( I } 

DO 521 J = 1 » MY 

JK= J+M I D— 1 

D( I , JK) =CTAU( I , J) 

521 CONTINUE 
GO TO 711 

718 DO 719 I =1 , MX 

MY = NU M ( I ) — M AD ( I ) + 1 
DO 921 J =1 , MY 

921 F(J)=CTAU( I,J) 

MY = MUM ( I ) 

L=MAD( I ) 

DO 920 J =2 , 1 00 
L = L- 1 

CT AU ( I , L ) = F ( J ) 

IF(L.EO.l) GO TO 922 
920 CONTINUE 

922 L=MAD ( I ) 

DO 923 J=L , MY 

923 CT AU ( I , J ) = D( I, J ) 

719 CONTINUE 
711 CONTINUE 

WR IT E ( 6 , 79 ) 

DO 14 1=1, MX 
WRITE (6, 80) I, MX 

WRITE! 6 ,81 ) DELSTR! I ), THETA! I ) , SHAPE ( I ) , EPSIL ! I ) 

WRITE! 6,82) 

MY = NUM ( I ) 

DO 13 J=1 , MY 
IFU.EQ.30) WRITE (6, 94) 

I F ( J • EQ . 60 ) WRITE! 6 , 94 ) 

I F ( J . E Q. 90 ) WR IT E (6 ,94 ) 

IF(J.EQ.120) WRITE (6,94) 

WRITE (6, 83 ) J, U( I » J ) ,URAT 10! I,J),P(I,J),CTAU!I,J) 

13 CONTINUE 

14 CONTINUE 
WRITE(6,84 ) 

WRITE (6, 85) 

DO 15 1 = 1, MX 

XD = 12 . O^'XPOS ( I J/DIAMTR 

WRITE! 6 , 86) I , XD , D E LST R ( I ) , D ELR AT ( I ) , TH ET A ( I ) , THR AT ( I ) 
1 , SHAPE! I ) , EPSIL ( I ) 

15 CONTINUE 

CALCULATE THE TURBULENT KINETIC ENERGY RESIDUAL AND MU 
DO 715 1=1, MX 
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MY=NUM ( I )-MAD { I ) + l 
MI D=MAD ( I ) 

DO 715 J = 1 » MY 
JK= J+MI D— 1 
CTAUCI , J)=D(I» JK) 

715 CONTINUE 

DO 230 1=1, MX 
J4 = 0 

M I D=M AD ( I ) 

MT = MUT ( I ) 

DO 231 J=MID,MT 
J 4= J 4+ 1 
MA= J-M I D+l 

231 F(MA)=(U(I,MT) t#2-U(l , J )**2 )*U I I , J ) 

CALL OTFE( Y, F, RAKE , J4) 

YF IT( I ) =RAKE ( J 4 ) *RHO 
230 CONTINUE 
MT=MUT ( 1 ) 

MI D=1 
J4 = 0 
1 = 1 

DE LX=X POS ( 1 + 1) -XPOS (I ) 

MY=NUM ( I ) 

M I D=M AD ( I ) 

DO 236 J=MID,MY 
MA=J-MI D+l 

F(MA)=CTAU( I,MA)*2.*( IUI I , J-l)-U( I, J + i) ) / I 2.*Y) ) 

236 J4=J4+ 1 

CALL OTF E ( Y,F, 3 AKE , J4) 

YFOTI = RAKEIJ4)-(YFIT I 1 J-YFITI I + 1))/DELX 
YFATI =RAKE I J4) + ( YFIT (1 )-YFIT I 1 + 1 ) J/DELX 
YFOT ( I ) =YF0T I / YF AT I 
DO 23 2 1=2, MX1 
MT=MUT( I J 

DELX=XPOS( 1+1 ) -XPOS ( I ) 

J4 =0 

MID = MAD( I ) 

DO 233 J=MID,^1Y 
M A = J - M I D+l 

F ( MA ) =C TAU( I » M A ) * 2 . * l ( U (I , J-l ) -U ( I , J + l ) ) / (2 .*Y ) ) 
233 J4= J4+ 1 

CALL QTFEIY,F, RAKE, J4) 

YFOTI =RAKE(J4)-(YFIT( I -1 J - YF I T ( 1+1) )/(2.*DELX) 
YFATI =RAKF(J4)+( YFITI I- 1 J - YF I T ( 1+ 1 ) ) / ( 2 . -DE LX) 
YFOTI I ) =YFOT I / YF AT I 

232 CONTINUE 
MT = MUT I MX ) 

J4 = 0 

I=MX 

MI D=MAD ( I ) 

DO 237 J = M I D , M Y 
MA= J— M I D+ 1 

F I MA ) = CT AU ( I, MA )*2 .*( I U I I , J-D-UI I , J+l) ) / ( 2.*Y) ) 

237 J4= J 4+ 1 

CALL OTFEI Y, p, RAKE , J4) 

YFOT I =RAKE I J 4 1-tYFIT ( I — 1 )— YFIT ( I ) ) /DEL X 
YFATI =RAKEIJ4)+I Y FI T I I -1 ) -Y F I T ( I ) ) / DELX 
YFOTI I ) =YFOTI/ YFATI 
I =1 

MT =MUT I I ) 

DELX = X POS I 1 + 1 ) -XPOS I I ) 

MI D = MAD I I ) 

J4 = 0 

DO 581 J=MID,MT 
MA= J— Ml D+l 
FI MA ) =U I I , J ) ** 2 
581 J4= J4+ 1 

CALL OTF E I Y , F, RAKE, J4) 

DER= (YFITI I + D-YF I TII )) /DELX 
AMU I I ) = D ER / RA< E I J 4 ) 

DO 2 34 I =2 , MX1 

DELX = XPOS( I + 1 ) -XPOS I I ) 
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MT=MUT ( I ) 

M I D = M A 0 { I ) 

J4= 0 

DO 23 5 J = M I D » M T 
MA=J-MID+1 
Ft MA } = U ( It J )**2 
235 J4= J4+ 1 

CALL QTFEt Y,F, RAKE , J4) 

DER= ( Y F IT ( I + 1)-YFIT ( I— 1) )/(2.*DELX) 

AMU (I ) =DER/RA<E( J4) 

234 CONTINUE 

WR I T E ( 6 » 96 ) 

DO 717 1=1, MX 

717 WR IT E ( 6 , 97 ) I » YFOT { I ) * AM U ( I ) 

CALCULATE EDDY VISCOSITY AND REYNOLD'S MIXING LENGTH 

WRITE (6, 87 ) 

DO 17 1=2, MX 

WR IT E ( 6 , 88 ) I 

MI N= ( (NUM( I) — 1 ) / 2 ) + 1 0 

MA X=NUM ( 11—10 

DO 16 J = MI N , MA X 

E 1 =ABS ( D(I,J)*0(II) 

E2 = ABS( U ( I , J+l) -U( I , J-l) ) 

E3 = 2 q*h 

EDDY ( 1 i J)=( El/ RHG) / (E2/E3) 

AMI XL ( I , J ) = SOR. T( El/(RHO*(E2/E3)**2) ) 

WRITE (6, 89 ) J , EDDY ( I, J), AM I XL { I, J) 

16 CONTINUE 

17 CONTINUE 

1000 CONTINUE 
999 STOP 
END 






THIS PROGRAM WAS DESIGNED TO INTEGRATE BY THE USE OF 
SIMPSONS RULE. 

INPUT: 

CALL INTEG(F t M,DELX, RESULT) 

F- ARRAY OF VALUES OF THEINTEGRAND 
M- THE NUMBER OF POINTS 
DELX- THE STEP SIZE 

RESULT- THE RESULT OF THE INTEGRATION 



# £ £ £ £ ^ i: £ £ # # ^ ^ ^ # # # # ^ ^ ^ ^ # # # 3}- * # 



SUBROUTINE I NT EG ( F , K , DE LX , RESULT ) 
IMPLICIT REAL*4 (A-H,0-Z) 
DIMENSION F { 1 ) 

C INTEGRATION BY SIMPSONS RULE 

K 1=K- 1 
A=F ( 1 ) +F ( K ) 

B = 0. 0 

DO 1 1 = 2, Kl, 2 
SUM=4.0*F< I ) 

B = B+ SUM 
1 CONTINUE 
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C = 0.0 

DO 2 1=3. K, 2 
SUM=2.0*F( I ) 

C=C+SUM 
2 CONTINUE 

RESULT=DELX 5 MA + B+C) /3 . 0 

RETURN 

END 



£ % £ sjc # i* A * £ # * 5>V # jj: A sjc sjc aje * # * fr si * *; * # £ # * # £ ••}: # A # £ * * $ ;jc # # * >;< £ 



THIS SUBROUTINE WAS DESIGNED TO FILL AND ARRAY KITH 
VALUE EXISTING AT AN INDEX IN THE ARRAY. 

INPUT: CALL FILL(M.N.U) 

M- AN ARRAY OF INTERGER VALUES OF THE INDEX 
N- A FUNCTION VARIABLE 

U- A DOUBLE DIMENSIONED ARRAY OF VALUES, U(I,J) 
WHERE THE INDEX M IS CONTAINED IN U. 






SUBROUTINE F I LL ( NUM , MX , U ) 
IMPLICIT REAL* 4 (A-H.O-Z) 
DIMENSION NUM( 1) ,U{ 12, 1) 
DO 2 J = 1 , M X 
M=NUM ( J J 
DO 1 I = M , 1 7 5 

1 U(J« I)-U(J.M) 

2 CONTINUE 
RETURN 
END 









THIS SUBROUTINE WAS DESIGNED TO TAKE AN ARRAY AND 
DETERMINE ITS MID-POINT BY MEANS OF AN ARBITRAY 
FUNCTION, AND RETURN THE INDEX OF THAT POINT. 

INPUT: 

CALL PRF IN { A, N. M ) 

A- THE ARRAY OF VARIABLES 

N- THE LIMIT OF THE VALUES TO BE CHECKED 

M- THE INDEX OF THE MID-POINT 






SUBROUTINE PRE I N { A , N, M ) 

DIMENSION All) 

DETERMINE MIDPOINT OF MAKE DEFECT 

DUM=1000 . 

N 1 =N— 2 
DO 1 I = 3 , N 1 
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B=A( 1-2 )+A( 1-1 )+A( I )+A( 1+1 ) +A( 1+2 ) 
C = AM INI ( B* DUM) 

IF(C.GE.DUM) GO TO 1 

DUM=B 

M= I 

1 CONTINUE 
RETURN 
END 






THIS SUBROUTINE WAS DESIGNED TO TAKE AN ARRAY OF 
VALUES AND SPLIT THEM AT AN ARBITRARY POINT, AND 
PLACE THESE VALUES IN ANOTHER ARRAY. 

INPUT: 

CALL SPLIT (A,B,M,N,JK) 



A- THE INPUT ARRAY 
B- THE OUTPUT ARRAY 
M- THE LIMIT OF THE SPLIT 
N- THE START OF THE SPLIT 

JK- THE NUMBER OF VALUES BETWEEN M, AND N 






SUBROUTINE SP L I T ( A , B , M , N , J K ) 
DIMENSION A ( L ) ,B( 1) 
IF(JK.EO.l) GO TO 1 
DO 2 I =N , M 
B ( I-N+l) =A ( I ) 

2 CONTINUE 
JK=M— N+l 
GO TO 5 

1 MD=M-N- I 
DO 3 1=1, MD 
I J = N-I + I 
B ( I ) =A ( IJ ) 

IF II J. LE.l ) GO TO 4 

3 CONTINUE 

4 JK=MD 

5 RETURN 
END 
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THIS FUNCTION WAS DESIGNED TO SELECT THE LARGEST VAL 
ARRAY, AND THE LARGEST VALUE IN THE TOP HALF OF THAT A 

I NPUT : 



FUNCTION UEE(U»MY,UM, MM » I , MT ) 

U- THE ARRAY TO BE SEARCHED 
MY- THE EXTENT OF THE ARRAY TO BE 
SEARCHED, AN INDEX 
UM- THE MAXIMUM VALUE IN THE ARRAY 
MM- THE INDEX OF THAT VALUE 
I- THE INDEX OF THE FIRST DIMENSION 
OF U { I , J ) 

MT- THE INDEX OF THE LARGEST VALUE IN 
THE UPPER HALF OF THE ARRAY 
UEE- THE VALUE OF THE LARGEST VALUE IN 
THE TOP HALF OF THE ARRAY 



& i'% £ if >•' jV £ Hz £ if 5 (f 



FUNCTION UE E ( J , MY , UM, MM, I , MT ) 
DIMENSION U( 12 , 1 ) 

UM = 0 . 

M=MM 

DO 1 J = 1 , M Y 
UI=UM-U ( I , J ) 

IF(UI) 2,2,1 

2 UM=U ( I , J ) 

F M - J 

1 CONTINUE 
UE E = 0 . 

I FtMM.GT .M) GO TO 3 
DO 4 J = M , M Y 
U I =UEE-U { I , J ) 

IF(U1) 5,5,4 

5 UEE=U( I , J ) 

MT = J 

4 CONTINUE 
GO TO 6 

3 UEE=UM 
MT = M M 

6 RETURN 
END 



THE FOLLOWING IS A PRINT OUT OF THE RAW PITOT-STATIC 
DATA FOR THE STUDY 



//G0.FT06F001 DD S P AC E= ( CYL , ( 1 , 1 )) 
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TAT I ON 4 PITOT STATIC DATA 
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AXIAL STATION 10UT OF11STATIONS 



APPENDIX F 
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15.500 111.86 0.99567 0.0 
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***SUMMARY OF INTEGRATED WAKE PROPERTIES*** 
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